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Abstract 



This thesis is devoted to a study of quantum fields at finite temperature. First, I consider 
Dirac fermions and bosons moving in a plane with a homogeneous static magnetic field 
orthogonal to the plane. The effective action for the gauge field is derived by integrating 
out the matter field. The magnetization is calculated, and in the fermionic case it is 
demonstrated that the system exhibits de-Haas van Alphen oscillations at low tempera- 
tures and weak magnetic fields. I also briefly discuss the extension of the results to more 
general field configurations. 

Next, the breakdown of ordinary perturbation theory at high temperature is studied. 
I discuss the need for an effective expansion and the resummation program of Braaten 
and Pisarski in some detail. The formalism is applied to Yukawa theory, and the screening 
mass squared and the free energy is derived to two and three loop order, respectively. 

The main part of the present work is on effective field theories at finite tempera- 
ture. I discuss the concepts of dimensional reduction, modern renormalization theory, and 
renormalizable field theories ( "fundamental theories" ) versus non-renormalizable theories 
("effective theories"). 

Two methods for constructing effective three dimensional field theories are discussed. 
The first is based on the effective potential, and is applied to field theory with N charged 
U (N) symmetric scalar coupled to an Abelian gauge field. The effective theory obtained 
may be used to study phase transition non-perturbatively as a function of N. The second 
method is an effective field theory approach based on diagrammatical methods, recently 
developed by Braaten and Nieto. I apply the method to spinor and scalar QED, and the 
screening masses as well the free energies are obtained. 
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Chapter 1 



Particles in External Fields 



1.1 Introduction 

One of the oldest problems in nonrelativistic quantum mechanics is that of a charged 
particle in a constant magnetic field. This problem was solved in 1930 by Landau [IJ, and 
the energy levels are called Landau levels. 

More generally, particles in external fields have been studied extensively since the early 
fifties, when Schwinger || calculated the effective action for constant field strengths in 
QED using the proper time method. The study of matter under extreme conditions such 
as very strong electromagnetic fields is of interest in various systems, and the applications 
range from condensed matter to astrophysics [3,4]. 

In the case of a constant magnetic field there exists another and perhaps simpler 
method for obtaining the effective action of the gauge field || . Integrating out the fermion 
fields in the path integral gives rise to a functional determinant, that must be evaluated. 
In order to do so we exploit the fact that the propagator equals the derivative of the 
effective Lagrangian with respect to the mass in the fermionic case, and with respect to 
the squared of the mass in the bosonic case. This requires the knowledge of the propagator, 
which can be constructed explicitly, since we know the solutions to the Dirac equation or 
the Klein-Gordon equation in the case of a constant magnetic field. 

Moreover, this method immediately generalizes to finite temperature and nonvanish- 
ing chemical potential. Thus, it becomes easy to study fermions and bosons at finite 
temperature and density. 

So far the gauge field has been treated classically. However, one may of course consider 
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quantum fluctuations around the classical background field. With the propagators at 
hand, one would then compute the vacuum diagrams in the loop expansion in the usual 
way. At the one loop level this implies a contribution to the effective action from the 
photons which equals ^^F"- ^ n 3+ Id, an d ^ 3 7 [ T in 2 + Id. This is the usual contribution to 
the free energy from a free photon gas at temperature T. Beyond one loop the evaluation 
of the graphs becomes difficult. Ritus has carried out one of the very few existing two- loop 
calculations in 3 + Id QED at T = 0, but finite density [|]. 

Many of the phenomena that have been discovered in condensed matter physics over 
the last few decades are to a very good approximation two dimensional. The most impor- 
tant of these are the (Fractional) Quantum Hall effect and high T c superconductivity @. 

Quantum field theories in lower dimensions have therefore become of increasing inter- 
est in recent years. Both systems mentioned above have been modeled by anyons, which 
are particles or excitations that obey fractional statistics. Anyons can be described in 
terms of Chern-Simons field theories [7-9]. 



Some ten years ago Redlich [10] considered fermions in a plane moving in a constant 
electromagnetic field. Using Schwinger's proper time method |2j to obtain the effective 
action for the gauge field, he demonstrated that a Chern-Simons term is induced by 
radiative corrections. The Chern-Simons term is parity breaking and is gauge-invariant 
modulo surface terms. 

External electromagnetic fields may give rise to induced charges in the Dirac vacuum 
if the energy spectrum is asymmetric with respect to some arbitrarily chosen zero point. 
The vacuum charge comes about since the number of particles gets reduced (or increased) 
relative to the free case. Furthermore, induced currents may appear and are attributed to 
the drift of the induced charges. This only happens if the external field does not respect the 
translational symmetry of the system. These interesting phenomena have been examined 



in detail by Flekk0y and Leinaas |1]J in connection with magnetic vortices and their 



relevance to the Hall effect has been studied by Fumita and Shizuya [T 

In this chapter we re-examine the system considered by Redlich We shall re- 



strict ourselves to the case of a constant magnetic field, but we extend the analysis by 
including thermal effects and we shall mainly focus on the magnetization of the system. 
For completeness, we also consider bosons in a constant magnetic field. Our calculations 
resemble the treatment given by Elmfors et al. [[| of the corresponding system in 3 + Id. 
However, interesting differences occur, mainly connected with the asymmetry in the Dirac 
spectrum, and we shall comment upon them as we proceed. 



1.2 Fermions in a Constant Magnetic Field 
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1.2 Fermions in a Constant Magnetic Field 

We start our discussion of particles in external fields by considering fermions in two 
dimensions in a constant magnetic field. 



1.2.1 The Dirac Equation 

In this subsection we shall discuss some properties of the Dirac equation equation in 
2 + Id. We also solve it for the case of a constant magnetic field along the z-axis. The 
Dirac equation reads 

(i<y"d lt -e'fA li -m)il> = 0, (1.1) 
where the gamma matrices satisfy the Clifford algebra 

Wrf} = 2<r. (i.2) 

In 2+1 dimensions the fundamental representation of the Clifford algebra is given by 2 x 2 
matrices and these can be constructed from the Pauli matrices. They are 

'-(:;)• --(!o). '-(;-")■ <-> 

Furthermore, in 2 + Id there are two inequivivalent choices of the gamma matrices, which 
corresponds to 7^ — > —7^. From Eq. (Q we see that this extra degree of freedom may 
be absorbed in the sign of m. These choices correspond to "spin up" and "spin down", 
respectively fTTH . 



The angular momentum operator o is a pseudo vector in 2 + Id, implying that the 
Dirac equation written in terms of these matrices does not respect parity This is 

no longer the case if the Dirac equation is expressed in terms of 4 x 4 matrices. These 
matrices can be taken as the standard representation in 3 + Id 

MS?)- ™ 

where we simply drop 74. This representation is reducible, and reduces to the two in- 




equivalent fundamental representation mentioned above ]T3 



In the following we make the choice 7 = a 3 , 7 1 = — ia 2 and 7 2 = — ia 1 . In this 
chapter we use the real time formalism and the metric is diag(l, —1, —1). For particles 
in a constant magnetic field, there are two convenient choices of the vector potential. 
These are = (0,By/2, —Bx/2) and A^ = (0,0, —Bx), and are termed the symmetric 
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and asymmetric gauge, respectively. In the first case the Hamiltonian commutes with the 
angular momentum operator and the solutions are given by Laguerre polynomials. In the 
second case we have [H,p y ] = and the solutions are Hermite polynomials (see below). 
We have chosen the asymmetric gauge and the Dirac equation then takes the form 



l et~ m 

- ieBx + #■ 

ox ay 



i-iz - ieBx + 9 



dx 



ay 



-i m -m 



V(x,*) = 0. 



;i.s) 



Here ?/>(x, t) is a two component spinor. Since the Hamiltonian commutes with p y , we can 
write the wave functions as 



^ K (x, t) = exp(—iEt + iky) 



f K (x) 
9k(x) 



(1.6) 



where k denotes all quantum numbers necessary in order to completely characterize the 
solutions. Inserting this into Eq. (1.5) one obtains 



E — m — £ + 

—E — m 



where 



£± = — id x =F i(k — eBx). 
The equation for f K (x) is readily found from Eq. ( |1.7| ): 

E 2 -m 2 -^)f K (x) = 0. 



The eigenfunctions of £+£_, provided that eB > 0, are 

1 , k 



7T 



1 4 cxp 



— ix j eB 

2 V eB J 



V2e~B( 



x 



eB' 



(1.7) 



(1.9) 



;i.io) 



Here, H n (x) is the nth Hermite polynomial. Furthermore, I Uik (x) is normalized to unity 
and satisfies 



iV2eBnI n _i 



£+I nik (x) = i^2eB(n + l)I n+1:k (x). 
Combining eqs. (|1.9|) and ( |1.10|) yields 

W) = I n ,k{x), E 2 = m 2 + 2eBn. 
The function g K (x) satisfies 



E + m 



f K {x), 



;i.ii) 



(1.12) 



1.2 Dirac Equation 
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implying that 



9k{x) = -iV2eBnI n -t 7 k{x). 
The normalized eigenfunctions become 



;i.i3) 



^S( x ; *) = exp(=R£„t + iky)* 



1 E„ ± m 



2E„ 



;i.i4) 



where n = 0,1,2,..., E n = y/m 2 + 2eBn and tp^(^.,t) are positive and negative energy 
solutions, respectively. Note that ^>q^(x, t) = and that we have defined J_i i a ; (x) = 0. 
The spectrum is therefore asymmetric and this asymmetry is intimately related to the 
induced vacuum charge, as will be shown in subsection ??. In Fig. a) we have shown 
the spectrum for m > and in Fig. [Tl] b) for m < 0. 



E 

A 



E 

A 



E 3 



E 3 
E 2 

Ex 



E = m 



E = m 



E-x 



E-x 



E-2 



E-2 
E-3 



a , 



») 



Figure 1.1: The energy spectra of Dirac fermions in the presence of a constants magnetic 
field, a) m > and b) m < 0. 
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The field may now be expanded in the complete set of eigenmodes: 



n=0 



dk r 



(+)( 



,(-)( 



;i.i5) 



Quantization is carried out in the usual way by promoting the Fourier coefficients to 
operators satisfying 



{K,k, &n',fc'} — 8n,n'$k,k', {d n ,k, d n i k ,} — 8 n ,n'&k,k'i 

and all other anti-commutators being zero. 



;i.i6) 



1.2.2 The Fermion Propagator 

In the previous section we solved the Dirac equation and with the wave functions at hand, 
we can construct the propagator. In vacuum it is defined by 



iS F (x',x) = (0| T[^(x',t / )*(x,0] |0>, 
where T denotes time ordering. By use of the expansion ( [L.15 ) one finds 



;i.i7) 



n=0 



dk 



iS F (x', x) = £ / — [ 6(t' - t)^(x', 0<*(*, t) - 9(t - 0<iV, t'Wnii*, t) 



The step function has the following integral representation 

Oj(t'-t) 



0(f - t) = — / du. 

2ni J uj — ie 



uj — le 

After some purely algebraic manipulations, we obtain 

E n + m 



1.181 



;i.i9) 



S F (x',x)ab 



4^2 £ / dkdu ~^~~ exp [ - *) + - y) ] x 



n=0 ' 
1 



CJ 2 - + 2£ 



Sab(n,uj, k). 



;i.2o) 



Here S a b(n, uj, k) is the matrix 



( m + u)) I ntk (x')I ntk (x) -iy/2eBnI ntk {x')I n ^ k {x) 
-i\ / 2eBnI n _ 1:k (x')I njk (x) (m - u)I n _ ltk (x')I n _ ltk (x) 



1.21) 



At finite temperature and chemical potential we write the thermal propagator as (see 
Ref. § for details) 

(S F (x', x))j3 >ft = S F (x', x) + S F ,fl (x', x). (1.22) 



1.2 The Effective Action 
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The thermal part of the propagator is 



oo . 

is*v^) = -£ / 

n=0 



2^ 



/+(^i2(x / ,ov'22(x,t) - /F(£„)<iV,o^(x,f) 



where 

;(+) 



i 



/i-V) = i-/r ; (-w 



(+)/ 



exp/3(cu + //) + !' 



1.23) 
(1.24) 



exp/3(u; — //) + 1 
This may be rewritten as 



i f 

Sp(x',x) = — y2 dkduo expik(y' — y) expiuj(t' — t)fp(uj)5(uj 2 — E\ — ie)S a b(n,uj, k). 

27T n=0 J 



Here 



f F (u) = 6{uo)F F + \u) + 9(-u)Ff\u). 



(-)/ 



;i.25) 



;i.26) 



As noted in Ref . |J , one is not restricted to use equilibrium distributions in this approach. 
Single particle non-equilibrium distributions may be more appropriate if e.g. an electric 
field has driven the system out of equilibrium. 



1.2.3 The Effective Action 

The generating functional for fermionic Greens functions in an external magnetic field 
may be written as a path integral: 



£(77,77,1^) = / VipVipexp 



i I d A x ( --F^F^ + ip(ip- m)ip -fjip + ipn 



. (1.27) 



The functional integral describes the interaction of fermions with a classical electromag- 
netic field. It includes the effects of all virtual electron-positron pairs, but virtual photons 
are not present. Taking this into account at the one- loop simply amounts to including a 
temperature dependent, but field independent term in C e s. 



The fermion field can be integrated over since the functional integral is Gaussian: 

Z^^AJ = det[i(i]p-m)}exp[i J d 3 x [--F^F^+J d 3 yrj{x)S F {x,y)n(y)]]. (1.28) 
Taking the logarithm of Z(r],fj, A^) with vanishing sources gives the effective action 



S, 



eff 



d 3 x 



77 TPV-v 



i Tr log [ i(ip— m) 



1.29) 
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Note that we have written log det = Tr log by the use of a complete orthogonal basis and 
that the trace is over space-time as well as spinor indices. Differentiating Eq. Q1.29Q with 
respect to m yields 

dCi 

— — = itrSWz, x). (1.30) 
dm 

The trace is now over spinor indices only. By calculating the trace of the propagator and 
integrating this expression with respect to m thus yields the one-loop contribution to the 
effective action. This method has been previously applied by Elmfors et al. || in 3+1 
dimensions. 



The above equation may readily be generalized to finite temperature, where we sepa- 
rate the vacuum contribution in the effective action 



C = C + d + = C + C 
where Co is the tree level contribution, and 



cff 



dCcs 
dm 



i tr 



Sf(x, x) + S F ,fl (x, x) 



[121) 



;i.32) 



Using eqs. (|1.20|) and Q1.21]) a straightforward calculation gives for the vacuum contribu- 
tion 

dk du 



trSpix, x) 



oo 

^" n=0 



u 2 - El + IS 



m 



^ n=l 

ieB 00 



ieB 



n—l,k 



X 



En 



n,k 



X) - 



An 



. , m ieB 
Integrating this expression with respect to m yields 



eBm 



d = — V Vm 2 + 2eBn + , 
2tt ^ 4vr 



1.33) 



;i.34) 



The divergence may be sidestepped by using the integral representation of the gamma 



function |0J and subtract a constant to make C\ vanish for B = 0, 

1 r°° ds 

C\ = 3- / — exp(— m 2 s) [eBs coth(e-Bs) — 1]. 

87T2 JO ,S2 



;i.35) 



This result calls for a few comments. We have chosen a gauge, where Aq = 0. However, 
we could equally well have chosen A to be a nonzero constant. This would give rise to 
an additional term in the effective action 

-,2 



5£ x 



m e 
\m\ 47r 



A B. 



;i.36) 



1.2 The Effective Action 
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This is simply the gauge dependent Chern-Simons term, whose existence first was demon- 
strated by Redlich [|10] . 



In the following, we shall only consider m > 0, except for subsection ??. Similar re- 
sults for m < can, of course, be obtained by the same methods. 

The finite temperature part of the effective action is calculated analogously using the 
thermal part of the propagator (|1.23|) . 



C f>,n 



TcB 00 

Y,[ M 1 + exp -(3{E n -fi)] + log[ 1 + exp -f3{E n + jj) 



2vr n=1 
TeB 

+— — log [ 1 + exp —f3(m — /x) 



2tt 



;i.37) 



Letting B — > it can be shown that one obtains the pressure of a gas of noninteracting 
electrons and positrons: 



= - / EdE 

mT 2 r 



2tt 



log[ 1 + exp —f3(E - fi) ] + log[ 1 + exp -f3(E + //) 

2tt 



Li 2 (-Xe- pm ) + Li 2 f-A- 1 e- /^m, 



Li 3 (-Ae- /3m ) + Li 3 (-A- 1 e- /3mN 



:i.38) 



Here, A = is the fugacity and Li n (x) is the polylogarithmic function of order n: 



oo ^/c 



fc=l 



;i.39) 



In the following we restrict ourselves to the case \i > 0. Analogous results can be obtained 
for /i < 0. 



In the zero temperature limit of CP^ one gets 



27r n=0 



1.40) 



where the prime indicates that the sum is restricted to integers less than (/i 2 — m 2 ) /2eB. 



Similarly, one may derive the density 

dCP>» eB ^ r 1 



P 



djjL 



+ 



2tt ^Uxp/^-//) +1 
e£ 1 
2tt exp /3(m — //) + 1 



exp/3(E n + /i) + 1 



;i.4i) 
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At T = this reduces to 



eB 
P= 2^ 



Int 



m 



2eB 



+ 1 



fj, > m, 



;i.42) 



in accordance with the result of Zeitlin ||16|| . /From Eq. (|1.42 ) one immediately finds that 
the density as a function of chemical potential for fixed magnetic field is a step function. 



This is intimately related to the integer Hall effect as noted in Ref. p7| . In Fig. |1.2j we 
have plotted the density as a function of chemical potential for low temperatures (dashes 
line: T/m = 1/100, solid line: T/m = 1/1000). One observes that the sharp edges get 
smeared out as the temperature increases. 



2Trp/eB 20 




/j,/m 



Figure 1.2: The density in units of eB/2ir as a function of \i/m for T/m = 1/1000 (solid 
line) and for T/m = 1/100 (dotted line). eB/m 2 = 1. 



1.2.4 Magnetization and the de Haas-van Alphen Effect 

In this section we study the physical content of the effective action which was obtained 
in the previous section. In particular we investigate a few limits to check the consistency 
of our calculations. 
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The magnetization is defined by || 



M 



dC 



cff 



OB 



The vacuum contribution to the magnetization is obtained from Eq. (|1.35f) 



M 1 



, 3 , 

i,7T2 JO 



00 expf— m s) 
as k 



S2 



escoth(ei?s) — 



e 2 £?s 2 



sinh (eBs) 



;i.43) 



;i.44) 



For the thermal part of the magnetization we find 



Te 



J- c- r 

— ]T log[ 1 + exp -/3(£ n - /i) ] + log[ 1 + exp -(3{E n + y) 

/.IT L 



2vr n=1 
Te 

+— log [ 1 + exp -/3(m - //) ] 
e 2 B ~ n r 1 



1 



1.45) 



2tt ^ E n L e xp/5(E n -/i) + 1 exp/3(£ n + /i) + l 
Magnetization at zero temperature. In the zero temperature limit Eq. Q1.45Q reduces to 

eBn 



2tt 



E 

n=0 



H-E n - 



E„ 



1.46) 



where the sum again is restricted to integers less than (/i 2 — m 2 )/2eB. The thermal part 
of the magnetization at zero temperature changes abruptly, when (y 2 — m 2 )/2eB increases 
by unity. Thus, oscillates wildly, in particular is the limit B — > not well defined. 

The strong field limit (B — > oo) of M^ ,Al is found to be 



2tt 



(/i — m). 



;i.47) 



In the weak I?-field limit (eB <C /i 2 — m 2 <C m 2 ) the vacuum contribution becomes 

e 2^> ^oo exp(— m 2 s) e 2 B 
as i 



Mi 



12vr 3 / 2 Jo 



S2 



12ir\m\ 



'1.481 



This agrees with the results of Ref. [18]. In order to get the strong field limit ( eB 3> m 2 ) 
of the vacuum contribution, we scale out eB and take eB — > oo in the remainder. This 



gives 



£i oc {eB)i Mi oc e^v 7 !?. 



;i.49) 



Vacuum effects contribute to the magnetization proportional to the square root of B. 
This should be compared with the corresponding result in 3 + Id, where the magneti- 
zation goes like -Blog(^) M. The thermal part of the magnetization was found to be 



12 



Particles in External Fields 



= e(/i — m)/27r. Hence, the vacuum contribution dominates, exactly as in 3 + Id. 
For /i^mwe see that the thermal part of the magnetization is nonzero. From Eq. ( |1 .42 ) 
one obtains p = eB/27c, so the nonzero magnetization is a consequence of the fact that the 
density increases as the magnetic field increases (since all particles are in the ground state). 



Magnetization at finite temperature. In Fig. |0] we have displayed the total magneti- 
zation as a function of the external magnetic field for different values of the temperature 
(fj>/m = 3/2, T/m =1/150 solid line, 1/50 dashed line, 1/5 dotted line). Fig. O] is a 



magnification of Fig. 1.3 in the oscillatory region. 



The fermion gas exhibits the de Haas-van Alphen oscillations for small values of 
the magnetic field. These oscillations have been observed in many condensed matter 
systems ||, and they were first observed experimentally in 1930 ||19|| . It is a direct 
consequence of the Pauli exclusion principle and the discreteness of the spectrum. 



We also note that the magnetization approaches a nonzero value as B 
specifically, in Ref. it is demonstrated that the limit equals 



Te 
2^ 



In 







(l + e-Km+M,^ 



0. More 



;i.50) 



Some comments are in order. It is perhaps somewhat surprising that the magnetization 
is nonzero in this limit. One should, however, bear in mind that the sign of m uniquely 
determines the spin of the particles (and antiparticles), implying that the system under 
investigation consists entirely of either spin up or spin down particles. This is not the case 
in 3 + Id, where the representations characterized by the sign of m are equivalent. By 
summing over ±m, or equivalently, by using four component spinors, one finds a vanishing 
magnetization as B goes to zero, exactly as in 3 + 1 dimensions. 

Finally, we have displayed the modulus of the vacuum part as well as the thermal 
part of the magnetization over a rather broad interval of values of eB/m 2 in Fig. O . 
(T/m = 1/5 and p/m = 1.1). The thermal contribution saturates for values of the field 
where the vacuum contribution starts to dominate. The reason is that all particles are in 
the lowest Landau level for high values of B, and that the energy of this level is indepen- 
dent of the magnetic field. 



High Temperature Limit (T 2 ^> m 2 ^> eB,p = 0). The high temperature limit is rather 
trivial. From a physical point of view, one expects that Cr^ approaches the thermody- 
namic potential of a gas of noninteracting particles of mass m. Indeed, in this limit, one 
may recover Eq. ( |1.38| ) by treating continuous variable. 
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Figure 1.3: The magnetization in units of em as a function of B in units of m 2 /e for 
different values of temperature, fi/m = 3/2. 
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Figure 1.4: Magnification of the oscillatory region in Fig. 1.4. 



14 



Particles in External Fields 



0.10 



0.08 - 



0.06 



Mj em 



0.04- 



0.02 



0.00 




Figure 1.5: Thermal and vacuum contributions to the magnetization for a Fermi gas. 
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We would also like to describe the system in terms of constant charge density. At 
zero temperature it is not possible to invert Eq. fll.42 ) to write the chemical potential 
as a function of density, since the step function is not one-to-one. However, \i can be 
interpreted as the Fermi energy at T = (as long as the highest occupied Landau level 
is not completely filled), so one can immediately write down the chemical potential as a 
function of density: 



'm 2 + 2eB Int 



,2vrp, 
eB' 



1.51) 



We should also point out that at T > there is a one-to-one correspondence between 
density and chemical potential (Fig. |1.2| ), so one can invert Eq. ( |1.41| ) numerically. 



We have used Eq. (|1.51|) to make a plot of the thermal part of the magnetization as 
a function of magnetic field at constant density and at T = 1/100. The resulting curve is 



displayed in Fig. 1.6 



0.05 



M/em 




Figure 1.6: The magnetization in units of em as a function of B in units of m 2 /e at 
T = 1/100 and constant density p/m 2 = 1. 

The de Haas-van Alphen oscillations are seen to be present for low temperatures and 
weak magnetic fields. Furthermore, it is seen that the magnetization is zero for large 
magnetic fields. This can be understood from the following physical argument: For large 
B-fields all particles are in the lowest Landau level and ji = m (recall that the degeneracy 
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increases linearly with B). The energy of the single particle ground state is independent 
of the external field (Eq = m), so increasing B cannot lead to an increase in C 13 '^, when 
the charge density (and therefore the particle density) is held constant. Hence, the contri- 
bution to the magnetization from real thermal particles vanishes in the strong field limit. 



1.2.5 Induced Vacuum Charges and Currents 

In this subsection we calculate the vacuum expectation value of the induced charge and 
current densities. Such calculations have been carried out in other contexts, e.g. in 
connection with magnetic flux strings (see Ref. jll)). We shall employ the most com- 
monly used definition of the current operator which can be shown to measure the spectral 
asymmetry relative to the spectrum of free Dirac particles. 



* a (x),(-f*(x) a ) . (1.52) 



Using the complete set of eigenmodes as given by Eq. ( |1.14| ), a straightforward calculation 
yields 

<><*» = -i£nr < L53 > 



Eq. Q1.53 ) is simply the Chern-Simons relation. It has previously been obtained by e.g. 



Zeitlin [16] using the proper time method. This result has the following physical inter- 
pretation: As we turn the magnetic field on, an unpaired energy level E = m emerges (in 
the case m > 0). The number of positrons therefore gets reduced relative to the free case. 
This can be interpreted as the appearance of electrons and results in a negative charge 
density. For m < a similar argument applies. 

A corresponding calculation of ( j (x) ) reveals that the induced current vanishes. This 
result should come as no surprise due to translational symmetry of the system. A non- 
vanishing vacuum current would arise in the presence of an external electric field and is 
then attributed to the drift of the induced vacuum charge. 



1.2.6 Conductivity and the Integer Quantum Hall Effect 



Let us next consider the conductivity. According to Ref. |2T[ the expression for the 
components of the conductivity a. t j can be expressed in terms of the polarization tensor 



1.3 Bosons in a Constant Magnetic Field 
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n^fcck): 



.an«(o,k) 



'i.i 



;i.54) 



and follows from linear response theory. Moreover, by considering the functional derivative 
of the effective action with respect to one may deduce that |21 

dp 



n Oi (/c o ,k^0) = 
Combining the above equations, one may infer that 



a 



ere,- 



dp 



6ije dB' 



[1.55) 



;i.56) 



Thus, the conductivity is Hall like. Using Eq. ( |1.41| ) and including the contribution from 
the induced vacuum charge, which was calculated in the previous section, we obtain 



,2 oo 



47 + 2^±i 



e 3 B 



E 



n=l 



n r 



+ exp (3(E n - p) 1 + exp (3(E n + p) 



expf3(E n - p) 



expf3(E n + p) 



2ttT ^ E n L [1 + exp (3{E n - p) ] 2 
Letting T — > one finds 

a = 



[1 +expf3(E n + p) 



2n 1 + exp f3{m — p) 
(1.57) 



e 2 e\ 
— + — Int 

47T 27T 



2 eJ B 



:i.58i 



We thus see that the conductivity is a step function for T = 0. The system therefore 
contains the integer Quantum Hall effect. This was also noted by Zeitlin |i~6| . The 
generalization of Zeitlin's result to finite temperature is new. 



1.3 Bosons in a Constant Magnetic Field 

In this section we focus the attention on bosons in a constant magnetic field. We calculate 
the effective action and derive the magnetization. We point out the differences between 
the bosonic and fermionic results. Finally, we generalize to constant field strengths and 
study pair production in a purely electric field. 



1.3.1 The Klein-Gordon Equation 



Let us for the convenience of the reader briefly discuss the solutions to the Klein-Gordon 
equation in an external constant magnetic field. It reads 

[D^ + m 2 }(f)(x,t) =0, (1.59) 
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where = <9 M + ieA^ is the covariant derivative and the metric is diag(l, —1, —1). We 
have again chosen the asymmetric gauge = (0, 0, —Bx) and assume that the wave 
functions are in the form 

0(x,t) = e ~ iEt+iky f(x). (1.60) 



The differential equation for f(x) then becomes 

i-U(*) = (E 2 -m 2 + eB)f(x) 



(1.61) 



where ^± were defined in Eq. Q1.8D and the eigenf unctions of were defined in Eq. ( |ITToD 
The normalized eigenfunctions of the Klein-Gordon equation are 



,eB, 



7T 



H„ 



k 

V2^B(x-—) 
eB 



;i.62) 



with corresponding eigenvalues E n = ym 2 + (2n + l)eB. The Klein-Gordon field can 
now expanded in the complete set of solutions: 



1 00 r dk 
= irY. / F a n,fc</>n,fc(x,t) +b* k <p* k (x,t) 
Air n=Q J E n 



;i.63) 



Quantization is carried out as in the fermionic case by promoting the Fourier coefficients 
to operators. The only nonvanishing commutators are 

Kfc, a} n < k \ = 4:7iE n 5 ntn <5(k - k'), [b n>k , b ] n , y \ = A-KE n 5 n ^5(k - k'). (1.64) 



1.3.2 Boson Propagators and the Effective Action 

The generating functional for bosonic Greens functions in an external magnetic field may 
be written as a path integral in analogy with the fermionic case 

Z(J,J\A I1 ) = JV(j)V^exp % j ' d 3 x (-^F^F^ + tf(D^ + m 2 )0 + J^ + ^J 

(1.65) 

We integrate out the bosons in the functional integral and get a functional determinant: 



Z(J,J\AJ = det[*(ZVX + m 2 )]exp [i J d 3 x [ - -F^F**" + j d 3 yJ^(x)A F (x,y)J(y)} 

(1.66) 

Taking the logarithm of Z(J, , A^) with vanishing external sources gives the effective 
action 



S, 



eff 



d 3 x 



itrlog i{D^D^ + m 2 



;i.67) 
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where we have written logdet = Tr log by the use of a complete orthogonal basis. The 
first term is denoted Co and is the tree level contribution. For a constant magnetic field 
we have £q = — Differentiating Eq. (|1.67|) with respect to m 2 yields 



dm 2 



—i trA F (x, x). 



(1.68) 



The next step is then to construct the boson propagator which in vacuum is defined as 

iA F (x',x) = (0| r[$(x / ,t / ) $t ( x ^)] |0>. (1-69) 
Here, T denotes time ordering as usual. By use of the expansion ( |1.63Q one finds 



A F (x',x) 



An 2 



n=0 



dk r 



(1.70) 

After some algebraic manipulations and using the integral representation of the step 
function, we obtain 



-I oo 

E 



A F (x',x) ■■ 
The trace then becomes 



4vr 2 — J or 



dkduj 



El + is 



exp -t) + ik(y' - y)} I n {x)I n {x'). (1.71) 



trAi?(x, x) 



4tt 2 



E 

n=0 



dkdu 



Or 



E 2 + 



is 



;i.72) 



a B ™ I 



47T 



n=0 



-En 



Integration with respect to m 2 gives the effective action: 

Be 



Ci = -— E V™ 2 + (2n + l)eS 

^ n=0 



1.73) 



Employing the integral representation of the T-function |T5[, we find 



efi Z" 00 ds 



-, 2 / 2 

57T2 JO S2 



sinh(e-Bs) eSs 



;i.74) 



The above expression has been rendered finite by requiring that C\ = for 5 = 0. This 
result is in accordance with the leading term in the derivative expansion employed by 
Cangemi et al. 



At finite temperature and chemical potential we write the thermal propagator as 

(A F (x',x))^ = A F {x\x)+A f3 /{x\x). (1.75) 
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The thermal part of the propagator is 



1 



iASV,,) = --2; 

* n n=0 J rj, ~ 



dk r 



fUEn)<l>nA^KA X )-fB(En)(f>nA^tk( X ) ■ ( L76 ) 



Here f^\uj) and /g ; (u;) are the bosonic equilibrium distributions: 



r(-) 



exp/3(co> — /i) — 1 



exp + /i) — 1 



Eq. ( |1.68|) is easily generalized to finite temperature. Writing C = C + C\ 
C + £ e ff, we have 



dC 



cff 



dm 2 



-i tr 



Straightforward calculations give the thermal part of the effective action 
C f>,» = ___ £ [ log [ i _ exp _p( En _ ] + log [ i _ exp _^ (jEn + /i) 



2vr 



n=l 



(1.77) 

(1.78) 
(1.79) 



The limit B — > is easily taken, and we find 

27T 7m 

mT 2 



DEE 



log[ 1 - exp -/?(£ - //) ] + log[ 1 - exp -p(E + n) 



2tt 



Li 2 (Ae^ m ) + Li 2 (A" 1 e- /3r ' 



^3 

2^ 



Li 3 (Ae^ m ) + Li 3 (A~V 



;i.8o) 



This is the minus the free energy for a gas of bosons, as expected. 



The limit T — > is trivial in the bosonic case. There is no Fermi energy, and all the 
particles are in the ground state. Hence 



C 



Am 



0. 



:i.8F 



Recall that we work with the grand canonical ensemble, so the above result implies that 
the pressure of the Bose gas vanishes. 

The high temperature limit equals the pressure of the Bose gas with B = as in the 
fermionic case. 



1.3.3 Magnetization 



The vacuum part of the magnetization becomes 

1 f°° ds c 



Mi 



n 2. / 2 

87T2 JO S2 



sinh(eSs) 



e 2 Bs cosh(e-Bs) 
sinh 2 (e5s) 



;i.82) 
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The thermal part is 



-L c r 
Z7T n=l 

e 2 5 ^ 2n + 1 



E 



2tt 2E n 



+ 



-exp[^(E n -/i)] -1 exp[/?(£ n + /i)] - 1 
Taking the weak field limit (eB <ti m 2 ) of Eq. (|1.82|) yields 

e 2_g ^oo _ expf— m 2 s) e 2 B 



1.83) 



Mi 



24tt 3 / 2 



247r|m| 



1.84) 



This is one half of the fermionic result. In the strong field limit we find that the magneti- 

3 1 

zation in the vacuum sector goes like e^B^. We have computed the vacuum and thermal 
parts of the magnetization numerically for the neutral Bose gas (T/m = 1 and fi = 0). 
The result is presented in Fig. [L7]. Note that we have plotted the modulus of the mag- 
netization. We see that the thermal contribution to the magnetization has a minimum, 
so the susceptibility changes sign. This was also observed in the corresponding system in 



3 + Id by Elmfors et al. []2"2" |. The system thus changes from diamagnetic to paramagnetic 
behaviour. We also note that the thermal part magnetization goes to zero as B — ► oo as 
can be seen from Eq. ( |1.83| ). This is in contrast with the fermionic case and stems from 
the fact that the single-particle energies increases with the magnetic field. 




Figure 1.7: The vacuum and thermal contributions to the magnetization for a Bose gas. 
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1.4 Bosons in Constant Electromagnetic Fields 

In the previous section we have considered the effective action for bosons in the presence 
of a constant magnetic field. The formula can rather easily be generalized to the case of a 
constant electromagnetic field, without doing any actual calculations. In 3+1 dimensions 
one can construct two independent Lorentz invariant quantities of E and B: 

- \f, v F^ = 1 -{E 2 - B 2 ), *F tlv F^ = E-B, (1.85) 

where *F flu = \e lllV(X ^F OL ^ . In 2+1 dimensions there is only one such quantity, namely 
\{E 2 — B 2 ). As a consequence of Lorentz invariance, one can simply let B — > v ' B 2 — E 2 
in the effective action. Doing so, and expanding the formula in powers of \JB 2 — E 2 , one 
finds 

e 2 {E 2 -B 2 ) 7e\E 2 -B 2 ) 2 
Cl ~ 48nm + 384(W + " ( 86) 

The first term in the above expansion may be removed by redefining the gauge field. This 
series expansion demonstrate the nonlinear behaviour of the electromagnetic field, which 
is inherit in quantum optics, but is absent at the classical electromagnetism. It is the 
two-dimensional (bosonic) counterpart of the famous Euler-Heisenberg Lagrangian found 
as early as in 1936 ||| . 

Effective Lagrangians and effective field theories is a major part of the present thesis, 
and will be discussed at length in later chapters. 



1.4.1 Pair Production in an Electric Field 

In this section we calculate the effective action in the presence of a constant external 
electric field E. In the case of | E | > | B | and in particular for | B \ = the effective action 
has an imaginary part. The physical interpretation of this imaginary part is an instability 
of the vacuum. This may be seen by appealing to the definition of the effective action as 
a vacuum to vacuum amplitude: 



out | in ) = e 



(out | in) 



-2lmr 



;i.87) 



Thus, if the effective action possesses an imaginary part the right hand side of the above 
equation shows that the probability that the system remains in the vacuum is less than 
unity and that pair production may take place. Letting B — * iE in Eq. ( |1.74j ) we find 



eE 
8^§ 



ds 



— exp{—m s) 

S2 



1 



sin(e-Es) eEs 



X 
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C\ now has poles along the real axis at s n = mr/eE. The integration contour should now 
be considered to lie slightly above the real axis. The contribution to the imaginary part 
of the effective action comes entirely from the poles and using the usual prescription to 



handle them [24], one finds: 



n=l 



712 



r— Ll3 

8tt 3 ^ 



— exp(— m 2 Ti/eE) 



(1.89) 



Eq. (|1.89| ) is then the probability per unit volume that the vacuum decays. Our result is 
very similar to that obtained by Schwinger ||. 



1.5 Concluding Remarks 

In the previous sections we have obtained the effective action for fermions and bosons 
in the presence of a constant magnetic field. In the fermionic case the most interesting 
findings were the de-Haas van-Alphen oscillations and the nontrivial limit B — » of the 
magnetization. We did not explicitly calculate the susceptibility, but it is straightforward 
to do so. However, for bosons we noted that the susceptibility changed sign. In the 
fermionic case, the susceptibility has been calculated and analyzed by Haugset in Ref. [pOfl , 
and it exhibits interesting structures as the magnetization itself. 

Our treatment could be extended in various ways. First, one should consider doing 
a two-loop calculation to incorporate the effects of virtual photons as well. At finite 
temperature, this is not trivial. The problem is infrared divergences due to the fact 
that the zeroth component of the polarization tensor, n Atl/ (fco,k), is nonvanishing when 
k = in the limit k — > 0. A careful study of the polarization tensor together with some 
resummation approach would be valuable. 



Secondly, one could consider more general field configurations than a constant mag- 
netic field, as has been done by Elmfors and Skagerstam in Ref. [25[ . 
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Chapter 2 



Resummation and Effective 
Expansions 



2.1 Introduction 



In the previous chapter we have seen how quantum field theory can be used to compute 
the expectation values of physical quantities in relatively simple systems. We also saw 
that the introduction of finite temperature complicated matters and made it necessary 
to compare scales m, \feB and T. In this chapter we turn our attention to finite T 
calculations in more ambitious field theories, where calculations are far less trivial than 
those of chapter one. 

Quantum chromodynamics (QCD) is today widely believed to be the theory which 
correctly describes strong interactions . As long as the number of quarks is sufficiently 
low, this theory is asymptotically free. This means that the coupling constant decreases 
with the energy scale, and that perturbative calculations can be carried out at high 
momentum transfer. Moreover, at long distances QCD becomes a strongly interacting 
theory, and lattice QCD indicates that there is a linear potential between two quarks in 
the strong coupling limit of QCD J26|. This potential is responsible for confinement, which 
is the fact that one never observes free quarks or gluons. The physical states are all colour 
singlets, which are bound states of quarks. These states are the familiar hadrons such as 
pions, nucleons and kaons. Lattice QCD also suggests that hadronic matter undergoes a 
deconfinement phase transition at sufficiently high temperature or high density |27| . This 



state of matter is simply a plasma which consists of free quarks and gluons. The early 
universe may very well have provided such extreme conditions (high temperatures), and 
so the study of the quark-gluon plasma is important in understanding the early universe 
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[28,29]. Furthermore, a quark-gluon plasma may also be produced in future colliders in 
heavy-ion collisions, and these experiments then makes it possible to study the deconfined 
phase directly |30 . 



There are several physical quantities of interest in a QCD plasma. One of them is the 
plasmon which is a longitudinal collective excitation. The real part of the longitudinal 
part of the gluon propagator gives the plasmon mass. The plasmon mass is equal to 
the plasma frequency which is the lowest frequency in the medium. The imaginary part 
of the longitudinal part of the gluon propagator yields information about the decay or 
lifetime of these excitations |HJ]. Historically, this quantity was extremely important for 
the development of a consistent perturbative expansion for quantum field theories at finite 
temperature. 

The first calculations of the damping rate 7 were based on conventional perturbation 
theory. However, it was soon realized that the result for 7 in naive perturbation theory 
is dependent on the gauge fixing condition, while the plasmon mass is gauge fixing inde- 
pendent (see Ref. p5| and Refs. therein). Moreover, in some gauges it also turned out to 



be negative, which has been interpreted as a plasma instability. 

The situation was, of course, unsatisfactory, since the damping rate is a physical quan- 
tity and if computed correctly it must be gauge invariant. The puzzle was around during 
the eighties, and led people to consider new propagators, which by construction are gauge 
fixing independent, and derive the damping rate using linear response theory. However, 
the results were mutually disagreeing and it was certainly not easy to discriminate be- 
tween the various values of 7. The resolution of this problem was given by Pisarski, who 
discovered that the one- loop result is incomplete ||33|| . The leading order result receives 
contributions from all order in the loop expansion (The first example of this was provided 
by Gell-Mann and Bruckner in the late fifties in nonrelativistic QED |33 ]. In order to 



get the leading order contribution to the free energy, one had to sum an infinite series 
of graphs, which were called plasmon diagrams). In other words, the naive perturbative 
theory breaks down and must be replaced by an effective expansion in which loop cor- 
rections are suppressed by powers of g. Pisarski was able to isolate this infinite subset of 
diagrams, which gave the leading order result, and resum them into an effective expansion 
that includes all effects to leading order in g. 

These results is a part of the so-called resummation program mainly due to Braaten 
and Pisarski and is the topic of this chapter. The effective expansion involves effective 



propagators and in nonabelian gauge theories also effective vertices, and is mandatory to 
use at high temperature, in order to obtain complete results. Resummed perturbation the- 
ory restores the connection between the number of loops in the loop expansion and powers 
of the coupling constant. Moreover, it cures the problem of gauge fixing dependence of 
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quantities that should be independent. In particular, Braaten and Pisarski, demonstrated 
the gauge invariance and also the positivity of the gluon damping rate. This settled the 
controversy of the gauge dependence of the damping rate, and initiated intense studies of 
thermal field theory. There are many important contributions, and also improvements of 
the original approach, and we shall comment upon them as we move along. 

A major application of resummed perturbation theory in recent years has been the 
calculation of free energies and effective potentials. The effective potential is an important 
tool in the investigation of phase transitions at finite T, for theories where some symmetry 
has been spontaneously broken at T = 0. Since the calculation of the effective potential, 
which is a static quantity, normally is carried out in the imaginary time formalism, it 
turns out that one only needs effective propagators for the bosons. Fermions need not 
resummation, and it is also sufficient to use bare vertices. The most important phase 
transition is the electro- weak phase transition in the standard model or extensions thereof, 
which took place in the early universe |36|. Its possible role for the baryon asymmetry 
that we observe today is basically a question of the nature of the phase transition. In order 
to generate any baryon asymmetry, the universe must have been out of equilibrium, and 
so the phase transition must be of first order. The electro- weak phase transition has been 
studied independently by several groups using resummed perturbation theory. Fodor and 
Hebecker |37j have calculated the two-loop effective potential in the standard model using 



Landau gauge. Arnold and Espinosa |38j have also investigated this phase transition as 
well as the Abelian Higgs model, applying a simplified resummation approach in which 
one uses an effective propagator for the n = bosonic mode only. The Abelian Higgs 
model is of interest in its own right, since this is a model of a relativistic superconductor 
(See Ref. @). 

The literature on calculation of free energies of quantum field theories at high temper- 
ature (T well above T c ) is now vast, and we shall comment upon only a few selected papers. 
Frenkel, Saa and Taylor were the first to push calculations beyond two loop in resummed 
perturbation theory. They computed the free energy to order A 2 in 4 -theory fjOfl , and 



Parwani and Singh have extended this to order A 5//2 ||41|| . Arnold and Zhai have computed 



the free energy in QED and QCD to fourth order in the coupling constant [42|. The free 
energy of high temperature QED has also been computed independently by Coriano and 
Parwani jKJ. In their papers Arnold and Zhai develop the machinery to deal with com- 



plicated multi-loop sum-integrals analytically, and this represents significant progress in 
perturbative calculations. Zhai and Kastening [03] have since extended the computations 
through fifth order. 



The outline of the chapter is as follows. In the next section we discuss the breakdown 
of perturbation theory and the resummation program of Braaten and Pisarski. In the fol- 
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lowing sections we apply resummed perturbation theory to Yukawa theory, and calculate 
the screening mass squared and the pressure to two and three loop order, respectively. 

In the Feynman diagrams, a dashed line denotes a scalar field and a solid line denotes 
a fermion. Our notation and conventions are summarized in the beginning of Appendix 
A and B. 



2.2 The Breakdown of Perturbation Theory 

We shall first list some definitions, introduced by Braaten and Pisarski in Ref. [ 55 |: 

• By high temperature (or hot field theories), we mean T>m, where m is any zero 
temperature mass. 



• A momentum (k , k) is called soft when both k and k = |k| are of the order gT. 

• A momentum (fc , k) is termed hard when at least one of the components is of the 
order T. 



We shall now discuss the breakdown of perturbation theory at finite temperature. Con- 
sider massless A0 4 -theory, for which the Euclidean Lagrangian reads 

C-E = \{d,<pf + ^p\ (2.1) 



The one-loop contribution to the two-point function is depicted in Fig. |2.1| and is inde- 
pendent of the external momentum. It gives a contribution 

Si(*b,k) = ^^. (2.2) 

This sum-integral is defined in Appendix A, and in dimensional regularization it is finite. 
The renormalized inverse propagator at one-loop is then 

T?(k ,k)=kl + e + ^. (2.3) 



From this expression we see that the one-loop correction to the full propagator is 
of the same order as the bare propagator for soft external momenta. Thus, the above 
calculation reveals that naive perturbation theory breaks down for k ~ V\T, and that 
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Figure 2.1: One-loop correction to the two-point function in scalar theory. 

one must use some kind of effective expansion in which loop corrections are down by 
powers of the coupling constant. As a part of the resummation program, we define an 
effective propagator which is the inverse of Eq. (|2.3|) . 

We may also write the one- loop correction as ^-^xthe tree amplitude, where P 
is the external momentum. More generally, loop diagrams that can be written in this 
way (where P characterizes the external momenta) are termed hard thermal loops. By 
definition then, hard thermal loops (HTL) are equally important as tree diagrams for 
soft external momenta. Moreover, the hard thermal loops receive their main contribution 
from a small region in momentum space, where the momentum is hard. 

Are there other hard thermal loops in A0 4 -theory than the tadpole? Or in other words, 
do we need resummed vertices as well as a resummed propagator? The answer is no, and 
below we shall demonstrate that the four-point function receives a one-loop correction 
which only depends logarithmically on the temperature at soft momenta. Hence, it suffices 
to use the bare vertex in perturbative calculations. The one-loop correction to the four- 
point functions is depicted in Fig |2]^ and the expression is. 




(2.4) 



N 




Figure 2.2: One-loop correction to the four-point function in scalar theory. 



The first term in the integrand is now written as 
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and correspondingly for the second one. Here, E p = p and so on. One will then encounter 
terms such as 

(2.6) 



fp [%LO n - E p ] [i(u n -u)- E p _ k ] ' 
where ui = k . The frequency sum is carried out by rewriting it as a contour integral in the 
complex plane |35J . The results are then expressed in terms of Bose-Einstein distribution 
functions, which is 

The term in Eq. ( |2.5| ) is then replaced by 

1 + n(E p ) + n(E p+k ) 



iu) — E p — E p+k 
and similarly for the others. This gives 



{21 



3, 2 f d s p 1 



[1 + n{Ep) + n{E p+k )}[- 



1 



-A 

2 J (27r) 3 AE p E p+k L L " ' '" v p/ ' '" v p ^ KJll ik ~ E p - E p+k ik + E p + E p+k ! 
-HEp) - n{E p+k )][- 1 - 1 ]]■ (2.9) 

Zfc ~ -&p + -&p+fc ^0 + — -c^p+fc J 

Let us now discuss the various terms in the above equation. The term that is independent 
of the Bose-Einstein distribution function represents the diagram at T = 0. It is logarith- 
mically ultraviolet divergent, and after renormalization it goes like ln-^-, where /z is the 
renormalization scale. For the T-dependent terms it is convenient to distinguish between 
soft and hard loop momenta. When both the external and the loop momenta is soft, the 
contribution to the integral is down by factors of A, since the soft momentum is the only 



scale in the integral [55 



The contribution from hard loop momenta can be estimated as follows. We use the 
approximations 

ik ± {E p + E p „ k ) w ±2E P ik ± {E p - E p „ k ) ^ik ±k cosd, (2.10) 

n{E p ) + n{E p _ k ) « 2n{E p ) n{E p ) - n{E p _ k ) w -n{E p ){\ + n (E p ))-^-, (2.11) 

E p _ k « E p - k cos 9. (2.12) 

These approximations are straightforward to derive by Taylor expansions in k/p. Here, 9 
is the angle between p and k. Plugging this into the first term in Eq. (|2.9|) reads 

3 2 r d 3 p n{E p ) 

8 J {2n)3p 2 {p-kcos9)' 1 ' ' 

The angular integral decouples from the radial integral and one finds 

-A 2 / *7^n(^W). (2.14) 
8 J {2n) 2 k y p +kJ K pJ K J 
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Noting that the distribution function cuts off the integral at p ~ T, we see that this term 
contributes only logarithmically in T to the in the one-loop diagram. This is also the case 
for the remaining terms, and we can conclude that loop corrections are down by powers 
of the coupling. 

There is yet another way to see the breakdown of perturbation theory due to infrared 
divergences. Assume that we would like to compute the screening mass of the scalar 
field. Generally, it is given by the pole position of the propagator and to leading order 
this is simply m 2 = AT 2 /24, which follows from the above calculations. Beyond leading 
order it becomes more complicated. Naively, one would expect that the contribution at 



next-to-leading order goes like AT and is given by the two-loop graphs shown in Fig. 2.3 



Figure 2.3: The two-loop graphs for the two-point function. 
However, this is incorrect. The first two- loop diagram reads 

-Tt Q pk^ (2 - 15) 

The term in which q = is linearly divergent^] in the infrared and assuming an infrared 
cutoff order m, the two-loop goes like A 3 ^ 2 T 2 . This diagram is the first in an infinite series 
of diagrams which are increasingly infrared divergent. They are called ring diagrams or 
daisy diagrams, and are displayed in Fig. |2.4|. One can easily demonstrate that they all 



contribute at A 3//2 to the screening mass, and hence the screening mass gets contributions 
from all orders in perturbation theory, analogous to the damping rate discussed in the 
introduction. 

Due to the similarity in this series of diagrams, one is able to sum it and although 
every term is IR-divergent, the sum turns out to be convergent both in the infrared and 
in the ultraviolet. Taking the symmetry factors into account one finds that a diagram 
with m loops yields a contribution 

f r 1 im-lj 1 

~i^ x %pA twr- (216) 



lr rhe setting sun diagram has a logarithmic divergence in the infrared. Hence, after curing the infrared 
divergence it contributes first at order A 2 In A. 
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Figure 2.4: Ring diagrams in scalar theory. 

We now restrict ourselves to the n = mode in the sum over q , since the other terms 
are down by powers of the coupling. Summing over m produces 

AT r ~ i fz l r = AT r 

Here, m 2 is the thermal mass at leading order. This integral is listed in Appendix B, and 
one finds 

XTm 

~ -g— • 2.18 

The nonanalyticity in A shows that it is nonperturbative with respect to ordinary per- 
turbation theory or that we receive contributions from all orders in perturbation theory. 
Alternatively, the infrared divergences that we have encountered reflects the fact we need 
to use an improved propagator in the perturbative expansion. 

It is straightforward to demonstrate that the improved propagator defined by the 



inverse of Eq. (2.3), actually resums this infinite set of diagrams. The improved propagator 
is then given 

A(fc ,k)= \ (2.19) 

However, in order to avoid double counting of diagrams, we must also subtract a mass 
term in the Lagrangian and treat this term as an interaction. We then split the Lagrangian 
into a free piece and an interacting piece according to 

C = I(^) 2 + imV, (2.20) 

Ant = -^m 2 2 + A 4_ (221) 

We can now recalculate the self-energy and demonstrate that it is really a perturbative 
correction to m 2 = 4^-. The self-energy is now given by the usual one- loop contribution 
as well as the new vertex which are shown in figure |2.5| . We get 

Si(fc , k) = -m 2 + ^" '' 



2 fp P 2 + m 2 
9 Xr 1 AT r 1 Am 2 r' 1 
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The prime indicates that the n = mode has been left out from the sum. Since this 
contribution is set to zero in dimensional regularization, we can still use the expression 
for the sum-integral listed in Appendix A. The self-energy is rendered finite by adding 
the mass counterterm and using the appendices, one finds 

M 2 = m 2 + £i(0,0) (2.23) 
XT 2 mT 
~24~ ~~ ~8T' 

which is down by vA, as promised, and also reproduces the leading part of the sum of 
the ring diagrams, given by Eq. ( [2.18| ). 




Figure 2.5: Diagrams contributing in the recalculation of the screening mass in scalar 
theory. 

Let us continue our discussion of the resummation program of Braaten and Pisarski 
by considering more complicated theories. The simplest theory containing fermions is 
Yukawa theory, which has been studied by Thoma in Ref. [[5J. The Euclidean Lagrangian 
is 

C E = \(d,<P) 2 + + m + gWt- (2-24) 

It is straightforward to show that one also needs an effective fermion propagator in Yukawa 
theory, and that is a common feature of all theories involving fermions. As in the pure 
scalar case one may show that the one-loop correction to the ijjipcj) also has a logarithmic 
dependence on T and it is not necessary to resum the vertex. One can demonstrate that 
this is the case for all other n— point functions in Yukawa theory, and hence we may 
conclude that only propagators need to be resummed. Let us take a closer look at this. 
The one-loop fermionic self-energy is shown in Fig. |2.6| and we have 



Zffa*) = 9% ptfp _* K)2 (2.25) 

The sum over Matsubara frequencies are carried out as in the bosonic case. The main 
difference is that the result involves both Bose-Einstein and Fermi-Dirac distribution 
functions, which is 

n{E p ) = —L^. (2.26) 
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Figure 2.6: One-loop fermion self-energy correction in Yukawa theory. 



After the appropriate substitutions and noting that we may neglect K in comparison with 
P in the numerator in Eq. fl2.25|) , we arrive at 



E/(fc ,k) 



«7o 2 
Sir 2 



d 3 p 



1 



{ 2,f i E r J l + n ^- h ^, h -E,,-E 
' 1 + [n(E p ) +h(E p _ k )][ 



p—k 



ik + E p + Ep^ k 
1 

iko — E p + E p _ k 

[ l 

ik - E p - E p _ k ik Q 
+ [n(E p ) + n(E p _ k )}[- 



iko + E p — Ep_k 



li 2 



d 3 p 



Pi 



(2tt)3 4E p E, 



p-'-'p—k 



[1 + n{Ep) - h{Ep_ k )} 



1 



Ep + E p _ k 
1 



iko + E p — E p ^ k iko — E p + E p _ k ~ 



(2.27) 



Now, one may infer that it is only the terms which involve the sum of two distribution 
functions that contribute of order T 2 ; The first term which is independent of the distribu- 
tion function, represent the T = contribution to the fermion self-energy and is therefore 
linearly divergent. The others are non- leading in T. Using the approximations above for 
hard loop momenta (the same approximations are valid for Fermi-Dirac distributions), 
we find that the angular integral decouples from the radial integral. The calculations are 
then straightforward and the fermion self-energy takes the form 



E/(fco,k) 



(2.28) 



Here, we have introduced the fermion mass m'j and the Legendre function of the second 
kind, Qo(x): 



9^ 
16 



m. 



Qo(x) 



In 



rx+ 1 



x 



The effective inverse fermion propagator then takes the form 

A-^fco.k) = -itf+iX f {ko,k). 



(2.29) 



(2.30) 



At this point we would like to comment upon the fermion propagator, which is given 
by the inverse of Eq. ( |2.30| ). As first pointed out by Klimov and Weldon, the effective 
fermion propagator has two poles [47,48]. The first corresponds to eigenstates where 
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helicity equals chirality and is the usual mode, known from T = 0. The second pole 
corresponds to eigenstates where helicity is minus chirality. This mode is a collective 
excitation and is occasionally referred to as the plasmino. 

Using the methods we have discussed in this chapter, the reader may convince herself 
that there are no other hard thermal loops in Yukawa theory. In particular, the ^i/jcf) 
receives a one-loop correction which depends on the temperature, only through logarithms, 
exactly as the quartic vertex in 4 -theory j|B . 



In QED, the only hard thermal loops are in the amplitude between a pair of fermions 
(fermion self-energy) and between a pair of photons (polarization tensor or the photon 
self-energy). The effective photon propagator is rather involved due to its non-trivial 
momentum dependence, in contrast with the local mass term in pure scalar theory. In 
QCD, there are hard thermal loops in all multi-gluon amplitudes, and also in the ampli- 



tude between a pair of quarks and any number of gluons ||35|| . Hence, effective vertices 
are required, too. The hard thermal loops have a remarkable property, namely that of 
gauge fixing independence. Klimov and Weldon were the first to show this property for 
the self-energy in QCD [47,48], and these results were extended by Braaten and Pisarski 
to all hard thermal loops by explicit calculations fl3"5" |. Moreover, Kobes et al. have given 



a general field theoretic argument of this property ||32|| . After the discovery of the gauge 
fixing independence, Braaten and Pisarski were able to construct an effective Lagrangian 
that generates the hard thermal loops for all amplitudes, and this effective Lagrangian 



is gauge invariant |49|]. Taylor and Wong constructed independently another equivalent 



effective Lagrangian with the same properties |50| . An effective Lagrangian that generates 
the hard thermal loops in Yukawa theory naturally also exists, and it reads 

2 - fdn 3 , 2 , f dtt d 2 , , 
C cS = C + ml^j / ^-^j + m 2 s (P / ^— 0. 2.31 

Here, we have introduced the four- vector K = (—i, k). The integral represents the average 
over the sphere. The fact that this effective Lagrangian generates the effective propagators 
follows easily from doing the angular integrals. In particular, the last term in Eq. Q2.31| ) 
is reduced to a local mass term. 

Above, we have seen how the improved expansion screens the infrared singularities 
that appeared in bare perturbation theory. However, in some applications it turns out 
that the HTL action does not screen all IR-singularities. A well-known example of this 
is the calculation of the electric screening mass in QCD beyond leading order |5l |. At 



next-to-leading order mass-shell singularities arise due to unscreened magnetic modes. A 
similar problem arises in scalar electrodynamics in the calculation of the scalar screening 
mass beyond leading order [52,53]. Another example of the breakdown of the original 
approach, is in the calculation of the production rate of real soft photons in a quark-gluon 
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plasma. For this problem, Flechsig and Rebhan have invented an improved effective action 
which removes these singularities [54fl. It can also be written in a gauge invariant way 
and reduces to the conventional HTL action, where the latter is valid. 



2.3 A Simplified Resummation Scheme 

In the previous sections we have discussed the resummation program of Braaten and 
Pisarski and demonstrated that we need to use effective propagators and, in some cases, 
effective vertices too. However, in the calculation of static quantities such as free energies 
(or effective potentials) and screening masses there exists a simplified resummation scheme 



due to Arnold and Espinoza p8| . The point is that for calculating Greens functions 
with zero external frequency, this is most conveniently carried out in the imaginary time 
formalism, without analytic continuation to real energies. We also know that in the 
imaginary time formalism, the Matsubara frequencies act as masses. For n ^ bosonic 
modes and fermionic modes they provide an IR cutoff of order T. Hence, for these modes, 
thermal corrections are truly perturbative (down by a factor of g), and it should be 
sufficient to dress the zero modes. Although it seems at first sight that the distinction 
between light and heavy modes may complicate things, it actually simplifies calculations 
a lot. We shall apply this approach in the next section to compute the screening mass 
and the free energy to order to order A 2 , Xg 2 and g A in Yukawa theory. 

Finally, we would emphasize that this simplified approach can not be applied in the 
calculations of dynamical Greens functions with soft external frequencies, as demonstrated 
by Krammer et al. in Ref. [[55 1. For instance, they show that it predicts an incorrect 



value of the plasma frequency at next-to-leading order in scalar electrodynamics, and 
they identify the problem to be that of an ambiguity in the analytic continuation to real 
energies. 
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As we have seen in the pure scalar case we must rearrange our Lagrangian according to 

£ = \^<P? + \m 2 5 kfl + + #V + 9W><I> ~ \m 2 5 kfi . (2.32) 

Here, the mass parameter m 2 is simply the bosonic self-energy at one-loop at zero external 
momentum. The relevant graphs are the first two Feynman diagrams in Fig. EO, and they 
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give 



Si(0,0) 
Using Appendix A we find 



2Tp P 2 



9~2L Tr 
'{P} 



pA 



m 



— T 2 _|_ ^ 2 rp2 



24 



6 



The resummed bosonic propagator is then 

i - 4,o 



A(fc ,k) 



K 2 



+ 



5, 



fc,0 



k 2 + w? 



(2.33) 



(2.34) 



(2.35) 



The screening mass is given by the location of the pole of the propagator at spacelike 
momentum |5T|. At the one- loop level, M 2 is then given in terms of the infrared limit of 



the self-energy function of the scalar field: 

M 2 = m 2 + £i(0,0). 



(2.36) 



Here, S n (0,k) denotes the nth order contribution to S(0,k) in the resummed loop ex- 
pansion. At the two-loop level we must take into account the momentum dependence of 
the self-energy function S(0, k). The screening mass is then given by 



M 2 = m 2 + 



Ei(0,k) + E 2 (0,k) 



k=im 



(2.37) 



Consider the leading order contribution to S(0, k), which is depicted in Fig. |2T7| . The 
second diagram is momentum dependent. Since the fermionic loop momentum is always 
hard, we can expand £i(0, k) in powers of the (soft) external momentum: 



Si(0,k) 



1-5, 



p,0 



P 2 



+ 



S 



p,0 



P 



77T 



9 2 i Tr 

T{P} 



P 2 (P + Kf 



— m 



XmT 
8vr 



+ 2k2 9 2 f {p} ^ + 0{k'jT 2 



(2.38) 



The sum-integral above is divergent and the divergence is removed by the field strength 
renormalization counterterm. To leading order we have 



1 - 



Thus, one finds 



Si(0,k) 



XmT 2g 2 k 2 , , A 

+ T ^(21n^+2 7s + 41n2 



57T 



(4tt) 



4ttT 



(2.39) 



(2.40) 



Here, A is the renormalization scale introduced by dimensional regularization (see Ap- 
pendix A) 
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Figure 2.7: Leading order contributions to the screening mass in Yukawa theory. 

Let us next consider the two-loop diagrams from the scalar sector. The first is inde- 
pendent of the external momentum k, and reads 



A 2 . L • 1 >\,a, S Pt o irl- 'V n <\ 



[- 



(7,0 Vqfl 



A 2 , 1 A 2 T 3 

+ 



4 Tpq L P 2 p 2 + m 2 ^ Q 2 q 2 + m 2 i 4 Tpq P 2 Q 4 3847rm 

\ 2 T 2 \ 2 mTr> 1 . 
128^ + lrfF' < 241) 

Note that here and in the following the prime indicates that the n = mode is left out 
in the sum (This does not affect the value of the sum-integral since the integral for the 
zero-frequency mode is set to zero). For the present calculation, the last term in the above 
equation is not needed, since it is of higher order in the couplings, and it is consequently 
dropped in the following. 

The second two-loop graph, the sunset diagram, is also the most complicated one. 
The terms in the sum for which at least one Matsubara frequency is nonvanishing are 
IR-safe, so that m is not a relevant infrared cutoff to order A 2 . Thus, we may put m = 
here. The remaining part where po — % — is infrared divergent and so we must keep 
the mass m. Hence, to order A 2 we can write 

_ 1 ~ flpo.ofrto.o f \ (o AO) 

6 Tpq P 2 Q 2 6 J Pq (p 2 + m 2 )(q 2 + m 2 )[(p + q + k) 2 + m 2 ] ' l ' ' 

Using the methods of Appendix C, one can demonstrate that the first term above is 
zero in dimensional regularization, and we are left with the second term. This integral 
is dependent on the external momentum k, and in order to calculate the screening mass 
consistently we must compute it at k — im. In order to see that this in fact is necessary, 
one can perform an expansion in the external momentum k, and verify that all terms are 
equally important for soft k ~ \^XT ~ gT. 

To this order we must also consider the tadpole with a thermal counterterm insertion. 
This is calculated the same way and one finds a contribution 

Am 2 r 1 — S Pi o 5 P:0 1 2 „ XmT 



+ 



P 2 P 2 + m 2 



16tt 
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Figure 2.8: Two-loop scalar diagrams. 




Figure 2.9: One-loop graph with a thermal counterterm insertion. 



Let us now turn to the the two-loop diagrams which come entirely from the Yukawa 
interaction. These are depicted in Fig. |2.10| and are all infrared safe when the mass 
m is set to zero. This implies that the leading order contribution, 0(g A ) can be found 
using the unresummed propagator. Moreover, due to the IR- convergence, m is not a 
relevant infrared cutoff and one may expand in the external momentum k. Thus, to 
0{g i ) it sufficient to consider the diagrams at vanishing external momenta. The first two 
diagrams obviously contribute equally and yield 
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4vf 



Tr 
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p e Q 2 (p + Qy 



*9 



1 



{P}Q P A Q< 
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The next diagram is treated in a similar fashion: 



9 4 t Tr 
J{PQ} 



P A Q\P + Qf 



±9 



1 



(2.43) 



(2.44) 



T { PQ } P2Q2(p + Q)2- 

This diagram actually vanishes in dimensional regularization, and this is demonstrated in 
Appendix C. 

The only mixed two-loop diagrams is infrared divergent when the mass is set to zero, 
so we keep the resummed propagator and find 



Tr{- 



+ 



r 1 — 5 n n 5„ 



2 T{P}Q ^P 2 (P + Q) 2 L Q2 g2 +m 2 

After the usual tricks we have 

\g 2 y.Ji' 1 T j: 1 



}■ 



1 T 4 1 

Y t*f{P}Q P 2 Q A + 2^T{P} P 2 " 



+ 



(2.45) 



(2.46) 



40 



Resummation and Effective Expansions 






Figure 2.10: Two-loop self-energy diagrams from the Yukawa sector. 

Here, the ellipsis indicate higher order terms, which can be dropped in the present calcu- 
lations. 



Figure 2.11: Mixed two-loop diagram contributing to the screening mass. 

Finally, we include the one-loop diagrams with counterterm insertions. These are 
depicted in Fig. |2.12| . At leading order we may again use the bare propagator and they 
contribute, respectively 



(2.47) 



was given in Eq. Q2.3S ), while the other renormalization counterterm reads [56 
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(2.48) 



The renormalization of the vertices are carried out by the replacements 
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2.4 The Screening Mass to Two-loop Order 
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Figure 2.12: One-loop graphs with wave function counterterm insertion in Yukawa theory. 

Evaluating Ex(0, k) at k — im and collecting our two-loop contributions, one finally 
obtains for the screening mass 
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This result here is new and it is easy to check that our result is renormalization group 
invariant by using the RG-equations for the couplings A and g 2 
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(2.52) 



Moreover, setting g = our result is in accordance with the one obtained by Braaten 
and Nieto using the effective field theory approach that we shall discuss thoroughly in 
the next chapters [)58| . We have also checked that this approach yields the same result 
also for g ^ (as well as the original approach by dressing all modes). We should also 
mention that the contribution from the one-loop with a thermal counterterm has canceled 
the two-loop contributions which individually were of the form X 2 T 3 /m and Xg 2 T 3 /m. 

Before closing this section we would like to make a few remarks. Firstly, consider 
the double-bubble, the mixed two-loop graph and the tadpole with a mass insertion. 
Individually, thse diagrams contribute to lower order than A 2 , Xg 2 and g A , but these 
cancel in the sum. This is of course necessary in order for resummed perturbation theory 
to work, and the reason for this is that the particular combination of these diagrams is 
infrared finite. This implies that to order A 2 , Xg 2 and g A , one may use the bare propagator. 
Above, we use the resummed propagator and calculated the diagrams, one by one, in order 
to demonstrate this cancelation explicitly. In the next section we shall use this observation 
to simplify the calculation of some three-loop graphs contributing to the free energy. 



Secondly, imagine that we would compute subleading contributions to the screening 
mass coming from the diagrams in the Yukawa sector. It is then mandatory to use the 
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resummed propagator, and we now show how to extract such contributions. Let us for 
simplicity confine ourselves to the second diagram which now reads 
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(2.53) 



As previously explained, this diagram is IR-safe in the limit m —>■ 0. The leading term 
is given by Eq. ( |2.44j ), and the subleading term is then found by subtracting the leading 
part (In this case the leading part accidentally vanishes, but that is besides the point). 
After changing variables, one finds 



^sub = -4m V 



T{P}Q P 2 (P + Q) 2 q 2 (q 2 + m 2 



(2.54) 



This expression is infrared divergent when the mass is set to zero. Hence, the second 
integral picks up its main contribution when q is of order m. Since P is always hard, one 
may to leading order in the couplings neglect q in comparison with P. Hence the integrals 
decouple, and the leading part of J su b is 
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(2.55) 



The other diagram yields a contribution at this order which is twice as large. The ultravi- 
olet divergences are canceled by renormalization of the quartic vertex in the third term in 
Eq. (|2.50p ). This contribution to the screening mass has a simple interpretation in terms 
of bare perturbation theory. It corresponds to multiple scalar self-energy insertions on the 
bosonic line in these graphs. These diagrams can be summed ad infinitum, and the situa- 
tion is indicated in Fig. [2.13| . In the last chapter we consider QED, and similar diagrams 
are present there (namely insertions of the photon self-energy). It is rather amazing that 
the corresponding contribution can be obtained by an almost trivial one-loop calculation 
in three dimensions! 





Figure 2.13: Infinite string of diagrams with scalar self-energy insertions. 



2.5 Free Energy in Yukawa Theory to order A 2 , \g 2 and g 4 
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2.5 Free Energy in Yukawa Theory to order A 2 , Xg 2 
and g 4 

In this section we shall compute the free energy to order A 2 , Xg 2 and g 4 using the methods 
from last section. The one-loop contribution is given by the following expression 



1„ f 1 „J1 , „o 9tt 2 T 4 Tm 3 



lnP 2 + -T/ — ^-22 lnP 2 = -^^ (2.56) 
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Figure 2.14: One- loop contributions to the free energy in Yukawa theory. 
The scalar two-loop yields 
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Note that to order we calculate we need not renormalize the second and third term in the 
above equation. The theta-diagram reads 
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Note again, that renormalization of the vertex is only necessary for the leading terms 
above. We therefore set Z 2 = 1 for the remaining terms. 

The one-loop diagram with a thermal counterterm insertion is 
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Note that we may combine the terms from Eqs. ( |2.57| ) and ( p.58|) which go like XmT 3 and 
g 2 mT 3 to obtain Hence it cancels the contribution from the one- loop diagram with 
a mass insertion. 
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Figure 2.15: Two-loop vacuum graphs. 




Figure 2.16: One-loop diagram with a thermal counterterm. 



Let us move on to the higher order contributions. At the three-loop level there are 
three diagrams which are infrared divergent. However, we also have a one- loop diagram 
with two mass insertions and a two-loop graph with one thermal counterterm insertion. 
These may formally be combined to a single graph, where the shaded blob denotes the one- 
loop self-energy for the scalar field minus the thermal counterterm. This is schematically 



displayed in Fig. |2.17| . The point here is that this particular combination is infrared finite. 
So we may use the bare propagator here, if we are to extract the leading contribution, 
which goes like A 2 , \g 2 and g A . For these diagrams, we can then write 



(2.59) 



where 



AE(P) 



1 

Q Q 



— 2 -3' 



Tr 

I 

to> Q 2 



(P + Q) 2 Q 2 



— m 2 S, 



Pofi 



t^2~ ~k + 2 9 2 i 

Tq Q 2 Tw\ Q 2 Jr 



{Q} Q 2 (P + QY 



m 5. 



(2.60) 



In order to proceed, we must distinguish between the light and the heavy modes. We 
first consider the heavy modes, and using the definition of the mass m 2 we easily find 
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Figure 2.17: Definition of the shaded blob. 

~ gA tp{ Q K } Q 2 K*(P + Qy(P + K) 2 ~ m2 fp{ Q } P 2 Q 2 (P + Q) 2 ' (2 ' 61) 
For the po = mode one obtains 

" ^{qk} Q*K\P + QYP + K f ■ (2 ' 62) 



We see that the fourth term in Eq. ( |2.61| ) combines with the term in Eq. ( |2.62 ) to the 



fermionic basketball. Note also that the last term in Eq. (|2.61|) cancels the last term in 



Eq ( |2.58| ), which follows from the fact that the fermionic setting sun diagram vanishes. 




Figure 2.18: Infrared finite combination of diagrams. 



The next task is to consider the infrared-safe three-loop diagrams. These are depicted 

in the propagators. 



in Fig. |2.19| . To leading order in the couplings, we can again put m 
The first one is simply the bosonic basketball: 
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The others stem from the Yukawa sector, and the first one reads 
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Here, £/(P) is the fermionic self-energy function defined in Eq. ( |2.25|) . 



(2.64) 



The second graph yields 
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Figure 2.19: Infrared safe three- loop diagrams in Yukawa theory. 

Finally, we have to include the two-loop diagrams with wave function renormalization 
counterterm insertions. The mass in the propagator is dropped for reason that should 
now be well-known. The first diagram is the double bubble, which reads 

(l-J^U i 



(2.66) 
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The theta diagram with wave function renormalization counterterms are displayed in 



Fig. [2.20| . They contribute 
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The renormalization of the vertices are carried out by using the expressions for Z\ and Z 2 
given earlier. Alternatively, all ultraviolet divergences at three loops can be canceled by 
renormalizing the coupling constants through the substitution A — > Z\\ and g 2 
in the two-loop diagrams. Here 
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Figure 2.20: Two-loop vacuum graphs with wave function counterterm insertions in 
Yukawa theory. 



Putting our results together, using Appendix A, we finally obtain the free energy 
through order A 2 , Xg 2 and g 4 : 
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This is the main result of the chapter on resummation and has not appeared in the 
literature before. It is easily checked that our result is renormalization group invariant as 
usual. Moreover, it coincides to order A 2 with previous results when g = |4(J . 
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Chapter 3 

Effective Field Theory Approach I 



3.1 Introduction 



The ideas of effective field theory or low energy Lagrangians have a rather long history 
and dates back to the early work of Euler and Heisenberg in the thirties, where they 
constructed an effective Lagrangian for QED, which could be used to compute low energy 
photon-photon scattering P^fl . In recent years, the applications of effective field theory 



ideas in various branches of physics have exploded. Most applications have been to systems 
at zero temperature (Refs. [59-65] and Refs. therein), but there is an increasing number 
of papers devoted to the study of effective field theories at finite temperature [66-76]. It is 
the purpose of this section to introduce the basic ideas of effective field theory in a rather 
general setting, and use the Euler-Heisenberg as a concrete example. We will also briefly 
discuss some major applications of effective Lagrangians at T = 0, that have appeared in 
the literature in the last couple of years. The introduction to effective field theories at 
finite temperature is deferred to the next section. 

The improved understanding of effective Lagrangian and the modern developments 
in renormalization theory have also led to some nice introductory papers to the subject, 
and we recommend the articles by Kaplan J77| , by Manohar |78| and by Lepage |79| . 



Now, what are the ideas of effective field theory and when can they be applied? 
Assume that we have a field theory which contains light particles of mass ~ m, and heavy 
particles of mass ~ M, where m « M so that one can speak of a mass hierarchy. Consider 
a process, e.g. scattering, which is characterized by energies far below M, so that no real 
heavy particles can be produced. It is then reasonable to believe that there exists an 
effective field theory for the light fields, which yields identical predictions for physical 
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quantities as the full theory in the low energy domain. This is an example of a general 
idea that pervades all physics; The detailed dynamics at high energies is irrelevant for the 
understanding of low energy phenomena. Of course, the high energy fields do affect the 
low energy world, but their effects may be fully absorbed into the parameters of the low 
energy Lagrangian. 

It is then the purpose of effective field theory methods to construct this effective 
Lagrangian which reproduces the full theory in the low energy domain, without its full 
complexity. Effective field theory ideas can, loosely speaking, be applied to any physical 
system with two or more distinct energy scales. One takes advantage of the separation 
of scales in the problem and treats each scale separately. This streamlines calculations, 
since we do not mix them. 

The effective field theory program can conveniently be summarized in the following 



points 80 



• Identify the low energy fields (the particle content) from which the effective La- 
grangian is built. 

• Identify the symmetries which are present at low energies. 

• Write down the most general local effective Lagrangian, which consists of all terms 
that can be built from the low energy fields, consistent with the symmetries. 

• The effective field theory can reproduce the full theory to any desired accuracy 
in the low energy domain by including sufficiently many operators in the effective 
Lagrangian. Specify this accuracy. 

• Determine the coefficients in the effective Lagrangian by calculating physical quan- 
tities at low energies in the two theories and demand that they be the same. This 
procedure is called matching. 

The fourth point above is in some sense a generalization of the Appelquist-Carrazone 



theorem [SI]: Consider correlators or Greens functions in the full theory with only light 
fields on the external legs, which are characterized by momenta k <C M. The Appelquist- 
Carrazone decoupling theorem then says that the Greens functions of the full theory 
can be reproduced up to corrections of order k/M and m/M by a renormalizable field 
theory involving only the light fields pTJ. By adding more and more operators to the 
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effective Lagrangian and tuning the parameters, Greens functions can be reproduced to 
any specified accuracy. 

Schematically, the effective field theory can be written as 

4* = 4> + £^> (3.1) 

n 

where we explicitly have isolated the renormalizable part, Cq, of C e s, and O n are operators 
of dimension n. The expansion in powers oik/M and m/M is referred to as the low energy 
expansion. 

The form of the effective Lagrangian can be very well understood in terms of the 
Wilsonian approach to the renormalization group . The starting point is the Euclidean 
path integral representation of the generating functional of the Greens functions in the 
full theory. We exclude the high energy modes in the path integral by using a cutoff A. 
This should be chosen to be much larger than M. The effects of the scale M have two 
sources; high energy modes of the light fields and all modes of the heavy fields. One way 
to isolate the effects of the scale M is to introduce a new cutoff A', so that m<C A'C M, 
and to integrate over all modes larger than A' for the light fields and over all modes 
for the heavy fields. The latter integration means that we actually eliminate the heavy 
fields from the path integral (integrating out the heavy fields). One is then left with a 
theory for the small momentum or long distance modes of the light field, and this has 
infinitely many terms. In this process the coupling constants get modified by the high 
energy modes, which is nothing but a renormalization of the parameters. The physics on 
the scale M is now encoded in the coupling constants. The fact that it is a local field 
theory follows from the Heisenberg uncertainty principle. The modes with momentum of 
order M are highly virtual and can only propagate over a distance of the order 1/M. So, 
at the scale m, this looks local and the effects of these virtual states can be mimicked by 
local interactions [ |79] . 

Let us take QED as an example. In this case the low energy field is the photon field. 
The symmetries are Lorentz invariance, gauge invariance, charge conjugation symmetry, 
parity, and time reversal. The next task is to specify the precision of the low energy 
theory. Since we are interested in describing photon-photon scattering we must include 
interaction terms in the effective Lagrangian. Let us for simplicity confine ourselves to 
the first nontrivial order in the low energy expansion. The most general Lagrangian to 
order /c 4 /m 4 (note that the electron mass is denoted by m, which is the heavy scale), 
satisfying the above requirements is the famous Euler- Heisenberg Lagrangian p3| : 
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Here, a, b and c are dimensionless constants. The parameter a is determined by considering 
the one-loop correction to the photon propagator in QED to leading order in k 2 /m 2 . The 
parameters b and c are determined by calculating the scattering amplitude for photon- 
photon scattering at energies well below the electron mass m in full QED and in the 
effective theory, and require that they be the same. These parameters can be written as 
power series in a, and the first term in this expansion was found by Euler and Heisenberg 
in Ref. [P5[ . Recently, b and c have been determined to order a 3 by Reuter et al. using 
string inspired methods p3] (See also ref. pi). 



Figure 3.1: Matching two-point functions at low energy. 




Figure 3.2: Light-by-light scattering in full QED, which can be mimicked by a local 
interaction on the scale A<m. 

After the constants a, b and c have been determined, we should be able to use it 
in the study of photon-photon scattering at low energy. At the tree level this presents 
no problem, but according to the traditional view on nonrenormalizable field theories it 
cannot be used at the loop level. Lets us consider this in some detail, and explain why 
the old view is incorrect. The operators -^(F^F^) 2 and ^p(*F flv F t11 ') 2 give rise to 7-7 
scattering, and they are nonrenormalizable, since their coupling constants have negative 
mass dimension. At the tree level, these operators reproduce the scattering amplitude in 
full QED up to corrections of order fc 6 /m 6 by construction, and all is well. The one-loop 



correction to photon-photon scattering is depicted in Fig. |373| and was first computed by 
Halter in Ref. f^ft 

The amplitude is dimensionless, and each vertex gives a factor 1/m 4 . Hence, the inte- 
gral must have dimension eight. If we use dimensional regularization, power divergences 

1 Halter's computation is correct, but it is does not provide the complete result to order k a /m s . The 
amplitude in the tree approximation is proportional to b and c, and consistency requires that these 
constants be determined to order a 4 . 
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Figure 3.3: One-loop correction to 7-7 scattering in the Euler-Heisenberg Lagrangian. 



are set to zero, and logarithmic divergences show up as poles in e [85]. 



The only mass scale in the integral is the external momenta fcj, and one will encounter 
divergent terms with more than four powers of fej. In order to render the amplitude finite, 
we must renormalize and absorb the divergences in the counterterms. However, there is 
no operator in the Euler-Heisenberg Lagrangian, which has dimension five or more. So, to 
get rid of the infinities, we must add one or more operators to the low energy Lagrangian 
which have the right dimension. Thus, we must introduce one or more coupling constants 
that must be determined from QED. Or, if we did not know the underlying theory, we 
had to carry out experiments to determine them. This argument can be repeated for any 
operator of a given dimension, and to any order in the loop expansion. This implies that 
£ e ff must in principle contain infinitely many operators, and that we have to determine 
infinitely many coupling constants. This is, of course, an impossible task, and has led 
people to the conclusion that the effective Lagrangian is useless and without predictive 
power beyond the tree approximation. This conclusion is not quite correct for the following 
reason: In order to reproduce the values of physical quantities to some desired accuracy, 
there is only a limited number of interaction terms in C e g, which has to be retained in the 
low energy expansion. The others are simply too suppressed by the heavy scale. So, if 
we are satisfied with finite precision, we only have to know a limited number of coupling 
constants in the effective Lagrangian, and this is naturally possible to calculate. This is 
relevant, because experiments are always performed with a finite precision. Hence, the 
effective Lagrangian is as good as any other quantum field theory, and can be used in 
practical calculations. 

The reader might nevertheless wonder: what is the point of constructing a low energy 
theory of photons? After all we do know the underlying theory, and full QED is probably 
not more difficult to use than C e s. This may be so in the scattering example, but it 
does serve as an illustration of the general philosophy. Moreover, the Euler-Heisenberg 
Lagrangian has other applications. A recent example can be found in Ref. [pO] , where 
Kong and Ravndal calculate the lowest radiative correction to the energy density for an 
interaction photon gas at temperature T m. The relevant vacuum diagram being the 
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photon double-bubble. The idea is again that of two widely separated mass scales m and 
T, and the correction goes like T 8 /m 4 . 

There is another important point we wish to make at this stage. We have argued 
that renormalizability is no longer a requirement of a useful and consistent quantum field 
theory. Why, then, is the very succesfull standard model renormalizable? This reflects 
the fact that the new physics first enters at a scale well above the one tested in this 
generation of accelerators. Moreover, there is today a large activity in the search for 
physics beyond the standard model, and future precision tests may very well reveal the 
presence of nonrenormalizable interactions which are signs of new physics. The bottom 
line is that before we have the final theory of everything (string theory or whatever) which 
applies to any energy scale, every field theory should be viewed as an effective field theory 
valid and with predictive power at a certain scale. Furthermore, even if we did know the 
theory of everything, it is unlikely that it will be of any use in low energy physics e.g. 
condensed matter. Reductionism in physics is more a principal question than a practical 
one. Thus, effective field theory methods are likely to be with us in the future as one of 
the most important tools for practical calculations. 

The full theory may not always be used at low energies, the most prominent example 
of this is quantum chromodynamics. QCD is a strongly interacting theory at low energies, 
and perturbation theory is useless in this domain. Hence, effective field theory methods 
are mandatory to apply. The essential ingredients in the construction of the effective 
Lagrangian is again symmetries and particle content. At the lowest energies only pions 
are present, and the symmetry is (an approximate) chiral symmetry (in addition to the 
space-time symmetries etc). The pions are approximate Goldstone bosons and in the 
chiral limit, they are massless. Chiral symmetry puts severe restrictions on the possible 
terms in the Lagrangian. Although the nonlinear sigma model has been around for many 
years, it is only in the last decade or so, it has been applied beyond tree level pTI . 



We shall give yet a few examples of effective field theory, which have received much 
attention in recent years. These are nonrelativistic QED and nonrelativistic QCD, which 
are effective field theories that are applied to bound states and were formulated by Caswell 
and Lepage |62j. This is an elegant alternative to the traditional approaches such as 



the Bethe-Salpeter equation, which do not take advantage of the nonrelativistic nature 
of the problem. Instead, it mixes contributions from different scales, making explicit 
calculations unnecessarily difficult. Let us for simplicity consider NRQED, which has 
been applied by several authors the last few years [63,64]. The form of the Lagrangian 
is again uniquely determined by particle content and symmetries. The fields present in 
NRQED are the two-component electron field, the electromagnetic field, and other fermion 
fields necessary for the actual problem. The symmetries are Galilean invariance, gauge 
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invariance, time reversal symmetry and parity. The coupling constants in NRQED are 
determined by the requirement that it reproduces full QED at low energy. This can be 
obtained by calculating scattering amplitudes in the two theories e.g. at threshold and 
demand that they be the same. This implies that the effects from nonrelativistic momenta 
are taken care of by NRQED, and that the effects of relativistic momenta are encoded in 
the parameters of the theory ||62|| . After the determination of the parameters in NRQED, 
one uses time ordered perturbation theory with the usual Schrodinger wave functions as 
the unperturbed states. 

Finally, there is general relativity. The modern ideas of effective field theory and 
renormalization theory have recently been applied to general relativity by Donoghue in a 
series of papers ||65|| . Conventional wisdom says that it is impossible to construct a mean- 
ingful quantum theory of general relativity, since the Lagrangian is nonrenormalizable. 
This apparent incompatibility of gravitation and quantum mechanics was considered as 
one of the greatest problems in theoretical physics. This is no longer so. As long as 
we are well below the heavy scale, which in this case probably is the Planck scale, it is 
perfectly possible to quantize gravity and it is a completely consistent theory at present 
energies. The fact that classical general relativity is in accordance with measurements 
simply reflects that the non-leading terms in the Lagrangian are strongly suppressed by 
the heavy scale M p . 

Of course there exist systems where the effective field theory program does not apply. 
This is the case if the energy scales in a system under consideration are not widely 
separated. Consider e.g. a hydrogenic ion with a large nuclear charge Z. This system 
is relativistic, which means the momentum of the electron is not much smaller than 
its mass. NRQED is therefore not a particularly useful approach and the low energy 
expansion converges very poorly. 



3.2 Finite Temperature 



Let us now begin our discussion of effective field theories at finite temperature. 



In the imaginary time formalism (ITF) of quantum field theory, there exists a path integral 
representation of the partition function f69 |: 



Z = / D</>exp 



dr I d 3 xCi 



(3.3) 



Here, is a generic field, and Ce is the Euclidean Lagrangian which is obtained from the 
usual Lagrangian after Wick rotation, t — > —it. Moreover, the boundary conditions in 
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imaginary time are that bosonic fields are periodic in the time direction with period (3 
and that fermionic fields are antiperiodic with the same period. 

The (anti)periodicity implies that we can decompose the fields into Fourier components 
characterized by their Matsubara frequencies: 

$(x,t) = /H[0 o (x) + 5> n (x)e 2 ™^], (3.4) 
*(x,r) = [3-^Mx)e m{2n+lW • (3.5) 

n 

The n = bosonic mode is called a light or static mode, while the n^O modes as well as 
the fermionic modes are termed heavy or nonstatic. In the ITF, we can therefore associate 
a free propagator 

A „(^) = _L_ (3.6) 

with the nth Fourier mode. Hence, a quantum field theory at finite temperature may be 
viewed as an infinite tower of fields in three dimensions, where the Matsubara frequencies 
act as tree-level masses of order T for the heavy modes, while the light mode is actually 
massless. 

In the preceding chapter we have seen that the scalar field acquires a thermal mass 
of order gT at the one-loop level. For the heavy modes this represents a perturbative 
correction which is down by a power of the coupling, and these modes are still character- 
ized by a mass of order T. However, the light mode is no longer massless, but its mass is 
of order gT. Hence, we conclude that we have two widely separated mass scales at high 
temperature, which are T and gT. The Appelquist-Carrazone decoupling theorem then 
suggests that the heavy modes decouple on the scale gT, and that we are left with an 
effective Lagrangian of the static mode. The process of going from a full four dimensional 
theory to an effective three dimensional Lagrangian is called dimensional reduction. This 
is the key observation and the starting point for the construction of effective field theories 
at finite temperature. These ideas were first applied in the eighties, and the main contri- 
butions from this period can be found in the papers of Ginsparg fBSfl , Jourjine [JTH] , and 
Landsman ||71|| . The parameters in the effective Lagrangian were determined by consid- 



ering the one-loop corrections from the nonstatic modes to the static n-point functions. 
In other words, only nonstatic modes circulate around in the loops and one speaks about 
integrating out the heavy modes. Thus, the effects of the scale T is now encoded in the 
parameters of the £ e ff, while the low energy effects should be fully accounted for by the 
effective field theory. 



Now, it was realized by Landsman that the effective Lagrangian does not completely 
reproduce the underlying theory [f7I| . This was taken as an indication that dimensional 
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reduction only takes place approximately. However, this apparent failure of dimensional 
reduction reflects the fact that only renormalizable field theories were considered. If 
one exploits the effective field theory program fully, and allows for nonrenormalizable 
interactions, dimensional reduction does take place in accordance with expectations 

Let us discuss the effective Lagrangian in the case where the underlying theory is 
A0 4 -theory. In the Feynman graphs below, light modes are indicated by dotted lines, 
while heavy modes are denoted by solid lines. 

First, we must identify the symmetries. We have a Z 2 symmetry — > — 0, which fol- 
lows from the corresponding symmetry in the full theory. There is also a three dimensional 
rotational symmetry. Hence, we can write 



£ eff = \{d^f + \m\K)<& + + ^ + hml^l + + - (3-7) 



Of the operators we have listed above, only the last one is nonrenormalizable. The 
coupling constants generally depend on the ultraviolet cutoff or the renormalization scale 
A. This is also the case for the field, but normally we shall suppress this dependence 
for notational ease. In our calculations, we use dimensional regularization, which, by 
definition sets the power divergences to zero, and where the logarithmic divergences show 
up as poles in e p5|pj 

Now, the coupling constants in the £ e fr are not arbitrary but determined by our 
matching condition, namely the requirement that the static correlators r^(0, k) in the 
full theory are reproduced to some desired accuracy by the correlators in the effective 
theory at distances R 3> l/T. This matching requirement was first introduced by Braaten 
and Nieto ||58|| , and independently by Kajantie et al. J73J. Below, we shall comment on 
the connection between the matching of Greens function and the old way of integrating 
out the nonstatic modes. 

The parameters in the effective Lagrangian are determined by ordinary perturbation 
theory in A, or g 2 , if we consider gauge theories. In the underlying theory this naturally 
corresponds to the following partition of the Lagrangian into a free part and an interacting 
part 



In the effective theory m 2 (A), A3 (A) and hi(A) are all of order A, while other constants 
are of order A 2 or even higher (although the operator 2 V 2 2 first contributes to the 

2 The power divergences are unphysical in the sense that the depend upon the regulator, and so 
dimensional regularization is a particular convenient choice. 



C 



hit 





(3.9) 
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screening mass at order A 5 / 2 and to the free energy at order A 3 ). This implies that we 
split the Lagrangian according to 



(£ cfr )o = \{di^) 2 (3.10) 



2 

(A^int = ^ 2 (A)0 2 + ^^ + ... (3.H) 



This way of carrying out perturbative calculations to determine the coupling constants 
was first introduced by Braaten and Nieto in Ref. pS |, and they refer to it as "strict 



perturbation theory" (see also chapter four). Now, we know from our previous discussion 
that strict perturbation theory is afflicted with infrared divergences, which are due to the 
masslessness of the fields. These divergences become more and more severe in the loop 
expansion, but we can nevertheless use it as a device to determine the coupling constants 
in £ c g. Below, we shall demonstrate that identical infrared divergences appear in the 
perturbation expansion in the effective theory. Thus perturbation theory breaks down in 
exactly the same way in the two theories, and the infrared divergences in the matching 
equations cancel. The coupling constants in C e g are only sensitive to the scale T at which 
perturbation theory works fine. We are namely integrating out the scale T, since the heavy 
modes have masses of this order, and the parameters encode the physics at this scale. So 
it does not matter that we make some incorrect assumptions about the physics on the 
scale gT (we dot use resummation) . However, when we use the effective three-dimensional 
theory in real calculations, we must take the screening effects properly into account. This 
amounts to including the mass parameter in the free part of the Lagrangian, and treat 
the other operators as perturbations. 

Let us now carefully demonstrate how this approach works in practical calculations. 
We shall outline how one determines the mass parameter at the two-loop level. We must 
start by determining the tree level Lagrangian. This is carried out by substituting the 
expansion of the scalar field in Eq. ( |3.4| ) into the path integral, and integrating over r, 
using the orthogonality of the modes. The expression we obtain contains terms which are 
made exclusively up of the static modes, terms that contain only nonstatic modes, and 
products of light and heavy modes. We can then read off the coefficients of the operators 
in £ e ff by comparing it to the part of the Lagrangian in the full theory that contains only 
the zero- frequency modes [f72]| . The reader may convince herself that the only nonzero 
coefficient in the tree approximation, is the parameter in front of the quartic coupling; 
A 3 (A) = AT. 



At the one-loop level in the full theory, there is one contributing diagram, namely the 
tadpole. This is depicted in Fig. where we have explicitly separated the contributions 
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to the static two-point function from the static mode and the heavy modes. One finds 
Here, the prime indicates that the n = mode has been left out from the sum. 




Figure 3.4: The tadpole graph, whose contributions from light and heavy modes have 
been separated. 

In the effective theory, the only relevant operators are m 2 (A)0Q and As(A)0q. The 
corresponding contributions to the two-point function are shown in Fig. |3.5| . The corre- 
sponding expression is 



T^(k) = e + m\A) + ^ ± (3.13) 



2 J P p 2 ' 

Demanding that these expressions be the same, determines the mass parameter to order 
A: 

2/AN A 3 (A) r 1 XT r 1 Xr'l 

Exploiting the fact that A 3 (A) = AT at leading order, we see that the second term on the 
left hand side cancels the first term on the right hand side. (Incidentally, this term is set 
to zero in dimensional regularization, since there is no scale in the integral, but that is 
besides the point). Hence 

- 2 (A) = ^^- (3.14) 



From Eq. ( 3.14| ), we conclude that the mass parameter is determined by the effects of the 



nonstatic modes circulating in the loop. This is actually a general feature of our matching 
procedure at the one-loop level; the first quantum correction to a coupling constant (which 
is the leading term in the expansion if the corresponding operator is not present at the 
classical level), is determined by the effects of the heavy modes in the loop. The matching 
procedure at the one-loop order coincides, not unexpectedly, with the original approach 



to dimensional reduction and effective field theories at finite temperature [71 . 



At two-loop order, it becomes more complicated. Integrating out the heavy modes 
in the above meaning of the word, implies that one considers the effects of heavy modes 
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Figure 3.5: The one- loop diagrams in the effective theory appearing in the matching 
procedure. 



in the loops. This is problematic, since it has been demonstrated by Jakovac that it 
generally produces non-local operators that cannot be expanded in powers of k/T |ff6| . 
Thus, it is difficult to construct a local effective field theory. A similar problem appeared 
in the study of QED many years ago. Ovrut and Schnitzer investigated QED with both 
light and heavy fermions P6| . They wanted to construct a low-energy theory containing 



the photon and the light fermion. Beyond one loop, they realized that one had to include 
both light and heavy fermions in two-loop graphs in order to obtain an effective theory, 
which reproduced the full theory at low energy. The solution to the problem at finite 
temperature is similar; One must be careful and consider diagrams with both static and 
nonstatic modes on internal lines. 



This is illustrated in Fig [3.6|, where we have displayed the graphs that contribute to 
the scalar self-energy function at two loops. The two-point function receives contributions 
from diagrams with only light lines, only heavy lines, with both light and heavy particles. 

Firstly, we would like to point out that the momentum dependence of the setting 
sun diagram is irrelevant for the present calculation. The momentum dependence of loop 
diagrams contributing to r^(0, k) gives rise to the renormalization of the fields in the 
effective theory. However, since this occurs at the two-loop level, this redefinition or 
renormalization first comes into play at the three- loop level (order A 3 ). 

We note that the two diagrams which contain only light modes (diagrams one and 
five) cancel against the corresponding graphs in the effective theory (see Fig |3.7|), exactly 
as at one loop. Moreover, the first of these diagrams in linearly infrared divergent, while 
the second has a logarithmic divergences in both the infrared and the ultraviolet. 
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Figure 3.6: The two-loop graphs for the two-point function, where the contributions from 
the light and heavy particles have been separated explicitly. 



The second diagram is zero, since the upper loop vanishes, while the third diagram is 
infrared divergent and reads 

A2T ^ 1 1 (3.15) 



4 J P p*r Q Qi 

However, this diagram is canceled by the diagram in the effective theory with a mass 
insertion. This graphs reads 

m 2 (A)A 3 (A) r 1 



pP 



4 ' 



(3.16) 



Consistency in the matching procedure, requires that we use the parameters m 2 (A) and 
A 3 (A) at leading order in A, and the cancelation then follows. The fourth diagrams gives 
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Figure 3.7: The two-loop graphs for the two-point function, in the three-dimensional 
theory. 

a contribution 

X 2 r 1 



4 Tpq P 2 Q 

while the sum of the sum of the last two diagrams is 



I* (3 - i7 » 



6 Tpq P 2 Q 2 (P + Qf 

Using the methods of appendix C, one can demonstrate that the above sum-integral 
vanishes in dimensional regularization. We can then summarize our discussion in the 
following matching equation: 

m . (A) + ftn . = ^^_^_^. (3,9) 

Note that we have included a mass counterterm on the left hand, which is necessary to 
cancel the divergence which is still left after the renormalization of the coupling constant^. 
The mass counterterm and the renormalization constant for A are, respectively p8|, 



T 2 \ 2 i\ 

Sm 2 = I * Z x = l + (3.20) 

24 167rV 32vr 2 e V ; 

Using appendix A, one obtains the mass parameter to second order in A: 

9, x AT 2 r A . A C'(-l)ii 
m (A) = 1 + - In 7^ + 2 + 2-9 — f . 3.21 



3 This divergence is related to the logarithmic UV-divergence of the setting sun graph in three 
dimensions. 
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The mass parameter is not renormalization group invariant, and this is a consequence of 
the A-dependence of the two-point function at two loop in the effective theory. 

We close this section by two comments. We have seen that the infrared divergences 
that arise at the two-loop in the underlying theory match the infrared divergence in the 
effective theory, as promised. This is a general feature of the present approach, although 
we shall not explicitly demonstrate this cancelation in the following calculations. 



Secondly, we saw that the sum of the last two diagrams in Fig vanishes, but 



individually they do not. This fact shows the difference between the matching procedure 
and the old way of integrating out the heavy modes. As previously noted, the latter 
method generates non-local operators. More precisely, the effects of the diagram with 
both light and heavy modes are incorporated by introducing a momentum dependent 
four-point vertex |76|. 



3.3 Spontaneously Broken Gauge Theories 

In the previous section, we studied the determination of the mass parameter in pure 
scalar theory at two-loop order. At the one-loop level this was straightforward, but 
the complexity increased at two-loop order. This will be even more dramatic in more 
complicated theories, which involves gauge fields. There exists a nice alternative to the 
direct evaluation of the Feynman graphs, namely the use of the effective potential. The 
main advantage of this method is that there are normally fewer diagrams involved in the 
calculations, and that symmetry factors are easier to figure out. This approach is due 
to Kajantie et al. [74||, who in series of papers study the construction of effective three 



dimensional field theories [68,73,74]. 

These effective theories are then used in perturbative studies and lattice simulations of 
the phase transition in spontaneously broken gauge theories. This includes investigation of 
the electro- weak phase transition in the standard model [I71J] or supersymmetric extensions 
thereof by Bodeker et al. fl3"6] , as well as SU(5) by Rajantie USTfl , and the Abelian Higgs 



model by Karjalainen and Peisa ||88|| . 

The idea is that the effective potential is the generator of one-particle irreducible 
Greens functions at zero external momentum. Now, assume that we split the Higgs field 
into a background field and quantum field in the usual way, and compute the effective 
potential in the full theory to some order in the loop expansion. The coefficients of 
0g/2 and (fio/24 then give the unresummed two and four-point functions, respectively, at 
zero external momentum. We then carry out a corresponding calculation of the effective 
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potential in the effective theory, and the above mentioned coefficients yield the same 
correlators in the three dimensional theory. Now, consider matching at one-loop. The 
contribution to any correlator from the n = mode is canceled against a corresponding 
contribution in the effective theory. Hence, the contribution from the n ^ modes to the 
one-loop effective potential in the full theory give, up to possible field redefinitions, the 
one-loop correction to any scalar Greens function. 

The determination of the parameters beyond one-loop is more complicated, since we 
must be careful with different mode contributions to the correlators. We shall make these 
ideas more precise in the next section. 

In the previous section we have seen that there are two mass scales in 4 theory, 
namely T and gT. This remark also applies to any Abelian gauge theory, where the scale 
eT corresponds to the scale of electric screening. In nonabelian gauge theories there is 
a third scale g 2 T, which is the scale of magnetic screening, or the inverse confinement 
radius. When the temperature is close to the critical temperature scalar fields have a 
mass of order g 2 T. This implies that we are faced with three scales in spontaneously 
broken gauge theories close to a phase transition, even in the Abelian case. Hence, it is 
useful to introduce the following definitions due to Kajantie et al. ]74|: 



• Superheavy modes. These are modes with masses of order T. The bosonic modes 
with n ^ as well as the fermionic modes are superheavy. 

• Heavy modes. These are modes with masses of order gT, where g is the gauge 
coupling. The temporal components of the gauge fields acquire masses proportional 
to gT and so these fields are heavy. For temperatures much bigger than T c the 
scalar masses is of the order gT and these modes are then heavy. 

• Light modes. These are modes with mass of order g 2 T or less. Near a phase tran- 
sition the masses of the scalar particles go like g 2 T and these modes are light. The 
spatial components of the gauge field are massless and so these modes are also light. 



Since we now have three different momentum scales, it is convenient to construct a se- 
quence of two effective field theories, which are valid on successively longer distance scales. 
This is a fairly straightforward generalization of the preceding discussion, so we summa- 
rize it in the following recipe: 

step 1 

Write down the most general Lagrangian consistent with the symmetries of the system 
containing light and heavy fields. The parameters are tuned so that the static correlators 
of the light and heavy fields in the full four dimensional field theory are reproduced by 
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the corresponding correlators in the effective theory, to some desired accuracy at distance 
scales R^> 1/T. The effective Lagrangian is valid for momenta k up to order gT, and the 
coefficients encode the physics on the scale T, which is a typical momentum of a particle 
in the plasma. The parameters are called short- distance coefficients. 

step 2 

This is an effective field theory where C' eS includes all operators that can be constructed 
out of the light fields, and which satisfy the symmetries of the system. The parameters are 
determined by demanding that the correlators of the light fields in the two theories match 
(to some required accuracy) at long distance (R ^> 1/gT). These coefficients give the 
contribution to physical quantities from the scales T and gT. Generally, the coefficients 
of operators in the two effective Lagrangians which involve only light terms differ. This 
renormalization of the parameters is of course due to the fact that we have integrated 
out the heavy fields and this difference encode the physics on the scale gT. This effective 
Lagrangian is valid for momenta up to order g 2 T, and the parameters are termed middle- 
distance coefficients. 



3.4 The Two-loop Effective Potential 

Now, let us apply the ideas of the previous section to a model which consists of iV 
charged scalars coupled to an Abelian gauge field. N — 1 corresponds to the Abelian 
Higgs model, that has previously been studied by several authors [38,39,89,90]. So the 
results presented here are a generalization of results that already appear in the literature. 
This generalization is fairly straightforward, but nevertheless very interesting. The point 



is that previous work on this model by Arnold and Yaffe [89], and by Lawrie p0| , using 
the epsilon expansion, indicate that the nature of the the phase transition depends on 
N. More precisely, for N larger than some critical N c ~ 365.9, the RG-equations have 
a nontrivial fixed point in coupling constant space in d = 4 — e dimensions. Such fixed 
points are taken as evidence for a second order phase transition U69"| . Thus from general 



considerations, one expects a first order phase transition for N < N c and a second order 
phase transition for iV > N c . However, perturbation theory normally breaks down for 
temperatures close to the critical temperature. Thus lattice simulations may be the only 
reliable tool in the determination of the order of the phase transition and our effective 
three dimensional field theory is the starting point of a nonperturbative study of the 
^-dependence of the phase transition. 



In the Feynman graphs a dashed line denotes the Higgs field. The Goldstone fields 
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are indicated by heavy dots, the wiggly line corresponds to the photon, and the ghost is 
denoted by ordinary dotted lines. The Euclidean Lagrangian is 

C = \f, u F^ + (P„$)t(P M $) - z/ 2 $t$ + ^($t$) 2 + £gf + £ gh . (3.22) 

Here D M = + ieA^ is the covariant derivative and $t = <I>2, •••> ®n)- ® * s the 
corresponding column vector. 

We perform the calculations in Landau gauge. The propagators and the gauge fixing 
term are, respectively 

A „„ (fc0 , k) = A H (k , k ) = (3.23) 

A cs (fa,k)= * 2 , £ ., = -i-fS^) 2 , a^O. (3.24) 

A z + 77^2 6 2a 

After the shift in the Higgs field, the tree-level masses are 

ml = -v 2 + ^0o, m 2 = -z/ 2 + ^ 2 , m£ = eVo- (3.25) 
2 6 

The effective theory of the zero modes consists of iV charged scalars coupled to an Abelian 
gauge field in three dimensions, in analogy with the full theory. The timelike component 
of the gauge field acquires a thermal mass and behaves as real scalar field. This field also 
couples to itself. We can then write 

+^|(A)0t0A 2 + ^-At + C gi + C gh + SC. (3.26) 

Here, SC represent all higher order operators consistent with the symmetries. The three 
dimensional gauge coupling is denoted by eg (A). 



The one-loop diagrams are shown Fig KM and they read: 



Vi = ~C(mi) + - l)C(m 2 ) + i(d - l)C(m v ) - ~C(0). (3.27) 

Here, we have defined 

C(m) = £ ln(P 2 + m 2 ) = £ In P 2 - ^ + m 2 ^' ^ - ±m 4 ^ + 0(m 6 /T 2 ) . (3.28) 



The first term is field independent and hence irrelevant for the present calculation. 
The second term comes from the n = mode, while the remaining terms arise from the 
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> o 




Figure 3.8: One-loop graphs contributing to the effective potential. 

n 7^ modes. These terms then give directly (up to the renormalization of the fields) 
the one-loop corrections to the mass parameter and the quartic coupling. The two-loop 
graphs are displayed in Fig |3.9| . Let us focus on the figure-eight diagrams which are in 
the form 



D ss (m 1 ,m 2 ) 



pq (P 2 + m 2 )(Q 2 + m 2 2 ) 
1 

pq (p 2 + m 2 )(q 2 + m\ 



+ T 



1 



1 



T 
T 



p p 2 + m 2 . JQ Q 2 + m\ 
1 

pq (p 2 + m 2 )(q 2 + m 2 , 

J £_L • 

P p 2 + m 2 .TQ Q 2 



+ 



+ T 



P p 2 + m 2 JQ Q 2 + ml 
1 

pq {P 2 + ml){Q 2 + ml) 
1 Ji' 1 



P p 2 + m 2 JQ Q 2 
1 



PQ P 2 Q 4 



(3.29) 



where the ellipsis indicates field independent terms as well as terms of higher order in 
m\ and m\. The first term which arise when both Matsubara frequencies vanish, is re- 
produced by the effective theory, and hence it is canceled in the matching procedure. 
Furthermore, the second and third terms are also canceled. There are (at least) two ways 
of seeing this. In order to obtain the effective two-loop potential resummation is required. 
The above-mentioned terms will then be canceled by the thermal counterterms |74| . Using 
the Feynman graph approach directly, such terms are canceled by one-loop graphs with 
a mass insertion in the effective theory. This was explicitly demonstrated in in the sim- 
pler scalar theory in section |3.2| . The last term above, where both zero-frequency modes 
have been removed, is then the only part of Dssijnx^mi) which contributes to the mass 
parameter. 



Consider first the theta diagram with three Higgs particles: 

1 



H(m 1 ,m 1 ,m 1 ) = W 
Tf 



pq (P 2 + m 2 )(Q 2 + m 2 )[(P + Q) 2 + m 2 } ' 



(3.30) 
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This diagram is IR-divergent in the limit m — > when both Matsubara frequencies vanish. 
Moreover, to leading order in the masses one can safely put m to zero, when at least one 
Matsubara frequency is different from zero (the diagram is IR-safe and the m is not a 
relevant scale at leading order. This is completely analogous to the the discussion in the 
previous chapter on resummation). Thus, we may write 

if (mi, mi, mi) — T 2 j — - — \-r- — ■ — — ^T+rf ^7^77^ +Q(mi) . 

^pg (P + mf)(p 2 + mf)[(p + q) 2 + m{\ Tpq P 2 Q 2 (P + Q) 2 

(3.31) 

As previously noted, the second term above vanishes in dimensional regularization, and 
so we are left with 

H{m 1 ,m 1 ,m 1 )=T 2 f — =_ \— + 0(ml). (3.32) 

V 1 Jpg G» 2 + m 2 )G» 2 + m 2 )[(p + q) 2 + m 2 ] v 11 v ' 

It is clear that this term is canceled by a corresponding term in the effective theory, when 
we match the two-point function. Hence, H(mi,mi,mi) does not contribute to the mass 
parameter m 2 (A). 

The next theta diagram is the graph with two vector propagators and one scalar prop- 
agator. By purely algebraic manipulations it can be written in terms of H(m 1 ,m 2 ,m 3 ) 
and double-bubbles. It reads 



D svv (m v ,m v ,m 1 ) = Y) 



pq (P 2 + m 2 v ){Q 2 + m 2 v )[(P + Q) 2 + m 2 ] 
(d — 2)H(m V) my, mi) 



4 

m 



Artiy 

771 r 

H — h- \H(mv,mi, 0) — H(m v ,m v , mi) 
m y L 

+2H(m v ,m v ,m 1 ) — -H(mv, mi, 0) 
-j- [/J55(my, my) - 2D S5 (m y , 0) 



^-[if(mi,0,0) + #(my, my, mi) — 2H(m v , mi, 0) 



— [-D S s(mv, my) + D ss (mi, 0) - D ss (m v , mi) 



2m y 



-%K,0)]. (3.33) 

From the above expression, we infer that the field dependent contributions from the figure- 
eight terms cancel. Hence, this diagram does not contribute to the mass parameter at 
this stage. 
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The final two-loop diagram can also be written in terms of -Dss(^i) ^2) and H(m\,m 2 , m$) 
and it is found to be: 

n , r . « (2P + QU2P + QU5, u -q,q u /Q^ 
Ussvymv, mi, m v 



(m\ - ml) 2 



pq (P 2 + m 2 )(Q 2 + m 2 v )[{P + Q) 2 + m 2 ] 
Dss(m v , mi) + D S s(m v , m 2 ) - D ss (m 1} m 2 ) 
+ (my — 2m\ — 2m\)H(my , m 1; m 2 ) 

if (my, m 1 ,m 2 ) - ff(m b m 2 , 0) 



my 



my 
m 2 . 



Hi 



V 



D S s(rn 2 , 0) + D S s(m v , m x ) - D S s(m v , m 2 ) - D ss (mi, 0) 

Dssirni, 0) + D ss (m v , m 2 ) - D ss (m v , mi) - D ss (m 2 , 0) 

(3.34) 



In contrast with the preceding diagram, this graph contributes to the mass parameter. 
More precisely, the second line above yields a contribution proportional to the mass of 
the vector particle. 



In terms of the sum-integrals defined above, the two-loop effective potential is 
A r 



V2 



24 



3D ss ( mi , mi) + (AN - l)D ss (m 2 , m 2 ) + (AN - 2)D ss (m u m 2 ) 



+ -(d-l)e 2 [D ss (m 1 ,m v ) + (2N - l)D ss {m 2 , m v ) 



H(m l} m 2} m 2 ) 



+e 4 (f)lD S vv(mv,rn v ,m 1 ) - -e 2 D ssv (mi,m 2 ,m v )}. 



(3.35) 



Finally, we must consider the counterterm diagrams. These are shown in Fig |3.10| . At 
the one-loop level, there is either a mass counterterm insertion or an insertion of a wave 
function counterterm. The corresponding integrals are 



ct 



1^ 5mj + 5Z^P 2 l/ 9Ar _ v O. Sm( + 5Z^P 2 



2Tp P 2 



m 



P 2 + 



m% 



Sm 2 v + 5Z A P 2 



P 2 + 



m 



(3.36) 



The counterterms above are those of the four dimensional theory, and they include those 
generated by the shift in the Higgs field. The mass counterterms read 

8m\ „ Q 3e 4 



c 2 u 2 \ 2(iV + 4)A 2 + 108e 4 , 2 



32vr 2 e 
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The wave function renormalization counterterms will be listed in subsection 13.5.1. 
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Figure 3.10: One-loop counterterm diagrams for the effective potential. 



3.5 The Short-distance Coefficients 
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3.5 The Short-distance Coefficients 

In this section we summarize our results of the discussion in the previous paragraph. We 
first calculate the field renormalization constants. We then use these and the effective po- 
tential to obtain the scalar mass parameter to two-loop order and the scalar self-coupling 
at the one-loop level. The other parameters in the effective theory are also determined to 
one loop order, and are found by computing Feynman graphs. 



3.5.1 The Field Normalization Constants 



In this subsection we determine the short-distance coefficients which multiply the fields 
in the effective theory. At leading order we have the simple relations: 
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(3.38) 



Beyond leading order this relation breaks down. At next-to-leading order one may read 
off the correction from the momentum dependent part of the two-point correlator, which 
is proportional to k 2 for $ and A , and k 2 5ij — kikj for A;. The coefficients in front are 
denoted by S'(0), IT 00 (0) and IT(0), respectively. Thus 
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For the gauge field the expression is 

Expanding to order k 2 and integrating by parts, we find 



i + n'(o)] 1/2 . 

(3.39) 
(3.40) 
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IL,(0,k) = — - — (Sij - kikj/k 2 )^ —. 

Correspondingly, one finds for the scalar field 

S'(0) = -3e^i,. 

After wave function renormalization, using 
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we find S'(0), n' 00 (0) and II' (0). The relations between the fields in 3d and 4d are 
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(3.45) 
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(3.47) 



3(4tt) 2 

We would like to emphasize that the first of the above relations is gauge fixing dependent, 
but the parameters of the effective theory are gauge independent. Secondly, taking into 
account the running of the fields in the full theory, we find that the three dimensional 
fields are independent of the renormalization scale A. This is related to the fact that there 
is no wave function renormalization in 3d. 




Figure 3.11: One-loop diagrams relevant for field strength normalization. 



3.5.2 The Coupling Constants 



The matching requirement for the quartic coupling yields the equation 

A3(A) = AT[1-S'(0)] 2 + Tr( 4 i ) ;, 1 , 1 (0). 



(3.48) 



The first term takes into account the different normalization of the fields, and the second 
is the one-loop correction, which may be directly read off from the effective potential 
Furthermore, the prime indicates as usual that we neglect the zero-frequency mode in the 
one-loop graph, and 4>\ is the Higgs field. This yields 
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(3.49) 



24vr 2 v 4vrT ' ' 4vr 2 - 

The couplings e 2 E (A) and h 2 E (A) are computed by considering the relevant correlators. 



The corresponding Feynman graphs are shown in Fig. |3.12| . The matching relations read 
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Again, the first term on the right hand side takes care of the different normalization of 
the fields, while the second term is the one-loop correction. These are given by 
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(3.52) 
(3.53) 



We notice that the term involving S'(0) in Eq. ( |3.50|) cancels against the first term in 
Eq. ( |3.52| ), after we have integrated by parts. This reflects the Ward identity. In the 
corresponding expression for /i^(A), there is not a complete cancelation, but we are left 
with a finite contribution. The net results are 
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(3.54) 
(3.55) 



Here, we have renormalized the electric coupling in the usual way to render the expressions 
finite. The charge renormalization constant is 



Z e 2 



1 + 



Ne 2 
48ttV 



(3.56) 




Figure 3.12: One-loop diagrams needed for the calculating the couplings e 2 E (A) and h E (A). 

We shall also compute the coefficient in front of the quartic self- interacting term Aq. 
At leading this can be found by considering the one-loop contribution to the four-point 
function for timelike photons at zero external momenta. The matching condition reads 

A A (A) = Tr^ oAoAo (0). (3.57) 

Note that the short- distance coefficient multiplying the 3c? fields does not affect the cou- 
pling at this order in e. We obtain 

rSlWO) = -12Ae 4 ^ -L + 96Ae^ j| - 96Ae 4 ^ ||. (3.58) 
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Figure 3.13: One-loop, four-point function with external timelike photon lines. 
This particular combination is finite, just as the corresponding one in QED f7l]: 



Ne 4 T 

X A (A) = -4— (3-59) 

Furthermore, A j a(A) first runs at order e 6 . 

The coupling constants A3 (A), e 2 E (A), hi (A) and A^(A) are renormalization group 
invariant. This can easily be seen by using the renormalization group equations for the 
scalar and gauge couplings at leading order [j3<| : 

de 2 Ne 4 

4 = < 3 - 60 > 

dX 1 
1 — = — 
d\x 24:71 

Thus, we can trade the scale A for an arbitrary renormalization scale fi. 



^ = _L_[(jV + 4)A 2 -18e 2 A + 54e 4 ]. (3.61) 



3.5.3 The Mass Parameters 

The mass of the timelike component of the gauge field to one-loop order follows directly 
from Eq. Q31ID 

9 , A N Ne 2 T 2 
m|(A) = — g— . (3.62) 

Normally, it is sufficient to know m% at leading order in the couplings, as above. In the 
next chapter we determine it at next-to- leading order (for N = 1), since we are specifi- 
cally interested in the electric screening mass. Its generalization to arbitrary N is not to 
difficult. 

The scalar mass parameter is determined from the two-lop effective potential. Schemati- 
cally we write 

m 2 (A) = -u 2 [\ - S'(0))] + Ei(0)[l - E'(0)] + E 2 (0). (3.63) 



3.6 The Middle-distance Coefficients 
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Here, Si(0) is the one- loop contribution from the nonstatic modes, and £ 2 (0) is the term 
that survives from the two-loop graphs. After renormalization, we are still left with a pole 
in e. This divergence is again a reflection of the logarithmic divergence of the two-point 
function in three dimensions. We cancel it by adding a mass counterterm, which thereby 
is determined to be 
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(3.64) 



This gives the mass parameter to two-loop order 
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(3.65) 



Note that the scalar mass parameter is explicitly dependent on the scale A. This depen- 
dence is necessary to cancel the logarithmic ultraviolet divergences in the 3d theory. 



3.6 The Middle-distance Coefficients 

We can also apply the effective potential to the three dimensional effective theory to inte- 
grate out the temporal component of the gauge field. The only relevant graphs are those 
with at least one Aq on an internal line. Diagrams which contain only light particles yield 
the same contribution on both side of the matching equation for any n-point function. 
The corrections to the parameters, found in this section are independent of N . 

The second effective field theory is identical to the previous one, except that the tem- 
poral component of the gauge field is left out: 

C' cS = = ^F l3 F lJ + (V^)\V^)+m 2 (A)^^+^^^^ 2 + C gf + C gh + 5C'. 

The gauge coupling is now denoted by cm(A), while the other parameters as well as the 
fields are barred in order to distinguish them from those of the previous section. The 
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Figure 3.14: Relevant graphs for integrating over the real scalar field A . 

dashed line denotes Higgs field, while the real scalar field is indicated by a solid line and 
the Goldstone particles are denoted by dotted lines. After a shifting the Higgs field in the 
usual way, the masses are 

m?(A) = m 2 (A) + ^M, (3.66) 

ml(A) = m 2 (A) + ^M, (3.67) 

6 

m 2 Ao (A) = m|(A) + e|(A)4 (3.68) 
The relevant diagrams are found in Fig. |3.14| and the contributions are 



V 2 = -t/W + <) - e%T^l ¥ - mW - , <)|(p - q)2 - m;] 

If 1 If 1 

+ 2 C|T 4 (p2 + mi f) )(g2 + m l) + ^ " 2 )e|T i p(? (p2 + m^,)^ + m 2 )' (3 ' 69) 

The bar on V2 is a reminder that contributions from the light fields to the two-loop ef- 
fective potential have been omitted. We have also multiplied by a factor T so that the 
effective potential has dimension four. 



Expanding in powers of e%(A)(fil/m 2 E (A) gives the one-loop contribution to the scalar 
mass parameter as well as to the quartic coupling constant. Close to the phase transition 
scalar mass goes like e 4 T 2 , and hence power counting implies that one can ignore the 
figure-eight graphs. 



3.6 The Field Normalization Constants 
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3.6.1 The Field Normalization Constants 

At leading order in the coupling constants the fields in the two effective theories are related 

as 

0(A) = 0(A), A i (A)=A i (A). (3.70) 

Now, the only one-loop diagram contributing to the scalar two-point function, where Aq 
circulates in the loop, is independent of the external momentum (this is the tadpole). This 
implies that there is no renormalization of the field 0(A) at leading order in the couplings, 
and the relation Eq. ( |3.70| ) still holds. Similarly, since Aq and Ai do not interact at leading 
order in e 2 , there is no renormalization of A either. This result is in contrast with the 



nonabelian case [74|, since A and A4 couples directly through the covariant derivative. 



3.6.2 The Coupling Constants 

The calculation of the gauge coupling e 2 M (A) turns out to be particularly simple, since Ai 
and Aq do not interact at leading order in e 2 E (h). Hence, the matching condition becomes 
trivial and we have 

4f( A ) = 4(A)- (3.71) 

This result is in contrast with nonabelian theories where Ai and Aq interact directly via 
the covariant derivative acting on Aj, and gives rise to a renormalization of the gauge 
coupling. Furthermore, e 2 M (A) is obviously independent of A. 

According to the discussion above, the quartic coupling gets a contribution at the one-loop 
level: 

The quartic coupling is also independent of the scale A. 



3.6.3 The Mass parameter 

Since the scalar fields does not get renormalized at the one-loop level, the contribution 
to the scalar mass parameter is directly given by the coefficients of 2 /2 from the sum of 
the one and two-loop graphs discussed above: 



m 2 (A) = m 2 (A) - 
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(3.73) 
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Here, we have canceled the pole term in the second integral in Eq. (|3.69| ) by a mass 
counterterm, 

and neglected m\ in comparison with my, using power counting. Notice again, that 
the mass parameter depends on the scale A due to the logarithmic dependence of the 
propagator at two loops in the effective theory. 

This concludes our calculations of the parameters of £ eS in terms of the tempera- 
ture, T, A, e 2 and the renormalization scale A. Finally, we mention, just for the record, 
that we have also calculated the parameters m 2 (A) and A3 (A) by evaluating the appro- 
priate Feynman diagrams. We obtain the same results, and it is a valuable check of our 
computations. 



Chapter 4 

Effective Field Theory Approach II 



4.1 Introduction 



In the previous chapter we studied in great detail the construction of effective field theories 
for the n = bosonic mode by matching Greens functions. Instead of explicitly dividing 
the loop corrections to static correlators into contributions from light and heavy modes, 
there exists a convenient alternative due to Braaten and Nieto E3. This approach is 
perhaps the most clean and transparent way of doing effective field theory, and calculations 
are greatly simplified compared to resummation methods, since one treats a single scale at 
a time. This is, of course, a common property of every effective field theory approach, but 
the simplification here is that one does not distinguish between static and nonstatic modes. 
The procedure of constructing the effective Lagrangian proceeds, though, essentially along 
the lines presented in the previous chapter. 

The starting point is the identification of the fields and the symmetries in the three 
dimensional Lagrangian. The parameters are tuned as functions of the couplings in the 
full theory, so as to reproduce the static correlators at long distances R ^> 1/T in the 
usual way. If there are two scales gT and g 2 T in the three dimensional theory as in e.g. 
QCD, one proceeds to construct a second effective field theory. We shall discuss this 
approach in detail in the present chapter. 

This method was first applied to A0 4 -theory by Braaten and Nieto in Ref. f58|. They 



computed the free energy to order A 5 / 2 , which was first obtained by Parwani using resum- 
mation methods ETJ. They also calculated the screening mass to order A 2 . Combining 



their results with renormalization group methods, leading logarithms of the coupling were 
summed. The latter result was new, and an improvement of the classic A 3//2 -result of Dolan 
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and Jackiw Later, Braaten and Nieto computed the free energy in QCD, through 



order g 5 [ 75j , and confirmed the resummation results of Zhai and Kastening j44j (which in 
turn was an extension of the g 4 -result of Arnold and Zhai |42|). Moreover, this approach 
provides a solution to the long-standing infrared problem in nonabelian gauge theories at 
high temperature to which we shall return in section [4.5| [75 . 



We shall apply these ideas to spinor as well as scalar electrodynamics. In QED, we 
obtain both the free energy and the electric screening mass squared to order e 5 . Our 
calculations reproduce previous results obtained by Zhai and Kastening [I3|, and inde- 
pendently by Coriano and Parwani ]92| for the free energy, as well as the screening mass 



first computed by Blaizot et al [p2 |. These results were obtained using resummation. 



In SQED we also compute the same quantities. In this case the screening mass 
squared is calculated to order e 4 and Ae 2 , and the free energy is derived to order A 2 , 
Ae 2 and e 4 . The former result confirms the calculations of Blaizot et al. who applied 
although they did not include a quartic self-interaction term for the 



resummation 52 



scalar field. The latter result is new and thus represents the present calculational frontier. 



4.2 QED at High Temperature 

In this section we discuss QED at high temperature (T ^> m) and the construction of 
the three dimensional effective field theory. The Euclidean Lagrangian of massless QED 
reads 

£qed = -F^F^ + ^7 M (<9 M - icA^ + £ gf + £ gh . (4.1) 

In this chapter all calculations are carried out in the Feynman gauge, but we emphasize 
that the parameters in the effective theory are gauge fixing independent. The gauge fixing 
term is then 

C gi =\{d,A,)\ (4.2) 
and the ghost field decouples from the rest of the Lagrangian. 

We call the corresponding effective three dimensional field theory electrostatic electro- 
dynamics (EQED), in analogy with the definitions introduced by Braaten and Nieto in 
the case of QCD JFBJ. The first task is to identify the appropriate fields and the symme- 
tries in EQED. It consists of a real scalar field coupled to an Abelian gauge field in three 
dimensions. The fields can, as usual, be identified (up to normalizations) with the zero- 
frequency modes of the original fields. In particular, the real massive field is identified 
with the n = mode of the timelike component of the gauge field in the full theory. Note 
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also that there are no fermionic fields in EQED, since the fermions decouple for reasons 
that should be clear at this stage. 

Now, £eqed must be a gauge invariant function of the spatial fields A i: up to the 
usual gauge fixing terms. This symmetry follows from the corresponding symmetry in 
the full theory and the Ward-Takahashi identity in the high temperature limit |[7I |. Since 



QED is an Abelian gauge theory, there will be no magnetic mass [Q. Moreover, the 
timelike component of the gauge field, A , behaves as a real massive self interacting scalar 
field. The fact that Aq may develop a thermal mass is a simple consequence of the lack of 
Lorentz invariance at finite temperature. Moreover, there is a rotational symmetry and a 
discrete symmetry A — > —A . The effective Lagrangian then has the general form 

£eqed = + \^A f + \ml{K)Al + ^-At + C gi + C gh + SC. (4.3) 

Here, A is the scale introduced in dimensional regularization. Furthermore, 5C represents 
all local terms that can be constructed out of A\ and Aq, which respect the symmetries 
of the theory. This includes renormalizable terms, such as qe{A.)Aq, as well as non- 
renormalizable ones like hE{A.)(FijFij) 2 . 



In Eq. (fO|), we did not include the unit operator. The coefficient of the unit operator, 
which we denote /^(A), gives the contribution to the free energy from the momentum scale 
T. So if we are interested in calculating the pressure we must determine it, as we determine 
other coefficients in EQED. If we are not, it is left out. Generally, this coefficient (as well 
the other parameters in EQED) depends on the renormalization scale A, in order to cancel 
the A-dependence which will arise from the calculations in the effective field theory. By 
including /s(A) in £eqed, we have two equivalent ways of writing the partition function 
in QED in terms of its path integral representation. In the full theory we have 

Z = J VfjVrjVA^Vip'D'ipexp^- dr J d 3 xC\, (4.4) 

where 7] denotes the ghost field. The result using the effective three dimensional theory is 

Z = e - fE{K)v J VrjVrt VA { VA exp [ - / d 3 x C EQED ] . (4.5) 

There is another physical quantity for a hot plasma, in addition to the free energy, which 
is of great interest, and this is the electric screening mass. This quantity gives information 
about the screening of static electric fields at long distances. The potential between two 
static charges in the plasma is normally derived in linear response theory f4"5(|, and reads 



r d 3 k 1 

y <*> = ^W" g + iMo,k) - (4 ' 6) 
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Here, Ilpf^k) is the photon polarization tensor. In the limit R — > oo, the potential is 
dominated by the pole in photon propagator. At leading order this pole is given by the 
infrared limit of il o(0, k), and the potential is thus a modified Coulomb potential with 
an inverse screening length or electric screening mass n o(0, k — > 0). This has led one to 
define the electric screening mass as the infrared limit of the polarization tensor ||45||: 



mi = n 00 (0,k^0). (4.7) 

This definition cannot be the correct one, since, beyond leading order in the coupling, 
it is gauge-fixing dependent in nonabelian theories [51,94]. Although U^ u (k ,k) is a 
manifestly gauge-fixing independent quantity in Abelian theories, the infrared limit is not 



renormalization group invariant, and is so a useless definition even here [51 



The electric screening mass is correctly defined as the the position of the pole of the 
propagator at spacelike momentum |51||: 



A; 2 + n 00 (0,k) = 0, k 2 = -mj. (4.8) 

This definition is gauge fixing independent order by order in perturbation theory, which 
can be proved on an algebraic level We also note that the two definitions normally 
coincide at leading order in the coupling constant. The above definition can be extended 
to other theories, e.g. 4 -theory. The polarization tensor is then replaced by the self- 
energy function of the scalar field, and the screening mass then reflects the screening of 
static scalar fields in the plasma. 

However, it turns out that one cannot calculate perturbatively the screening mass 
beyond leading order in nonabelian gauge theories using Eq. ( |4.8j ) |5Tf . The problem is a 



linear mass-shell singularity. This signals the breakdown of perturbation theory, and calls 
for a gauge-fixing independent and nonperturbative definition of the electric screening 
mass P5[]. 



In Abelian gauge theories the above definition is equivalent to defining the Debye 
mass as the correlation length of equal-time electric field correlation function |95| 



(E(x) • E(0)) ~ e- msX /x 3 , x^oo. (4.9) 



Here x = |x|. Unfortunately, the definition Eq. ( f4.9|) is a poor one in nonabelian theories, 
since E is no longer a gauge invariant quantity. The above considerations have led Arnold 
and Yaffe to define the electric screening in terms of Polyakov loops p5| . We shall not 



pursue this any further, but stick to the definition based on the pole of the propagator. 

1 The pole position is also independent of field redefinitions. Since the relation between the fields in the 
underlying theory and the effective theory can be viewed as a field redefinition, and since the screening 
mass is a long-distance quantity, one can use the effective theory to compute it. 



4.3 The Short-distance Coefficients 
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4.3 The Short-distance Coefficients 

In this section we determine the parameters of EQED. In the Feynman graphs a solid line 
denotes fermions, the photon is a wiggly line, while the ghost is indicated by a dotted 
line. In the effective theory the same conventions apply to the gauge field and the ghost, 
while the real scalar is denoted by a dashed line. 

As we have discussed in some detail in the previous chapter, strict perturbation theory 
is ordinary perturbation theory in e 2 . In full QED this corresponds to the usual partition 
of the Lagrangian into a free and an interacting part 

(£ QED )o = |^i + #V> + % + £ gh , (4-10) 

(^QED)int = -le^JJilp (4.11) 

The effective Lagrangian is split the following way: 

(£eqed)o = ^-^i-^i + 2^ iJ ^ ^ ^gf + (4-12) 

(£EQED)int = \<(ml+^< + 6£. (4.13) 



4.3.1 The Coupling Constant 

We need the coefficient in front of the quartic term Aq. This coefficient can be found 
by considering the one-loop contribution to the four-point function for timelike photons 
at zero external momenta, and was first obtained by Landsman in Ref. |[71|| . Hence, the 
matching condition is 

Xe(A)=TT^ AoAoAo (0). (4.14) 

The corresponding Feynman graph is displayed in Fig. ^T]. It requires the calculation 
of the following sum-integral, which is infrared safe since only fermionic propagators are 
involved: 



P} L 

4 r r 32pg 32p 2 , I 



70p70p270p27°p2 



The sum of these three integrals is finite in dimensional regularization and so the net 
result turns out to be 

2p 4 T 

X e (A) = (4.16) 

7T Z 
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This coefficient is independent of the renormalization scale A to order e 4 . Notice the sign 
in front of it, which is the opposite as in SQED, as we saw in the previous chapter. In 
SU(N) coupled to fermions, the sign depends on the ratio N/Nj, where Nf is the number 
[nfl . However, it is not large enough to shift the minimum of the effective 



of flavours 



action to a non-zero value of An 



Figure 4.1: One- loop four point function with external timelike photon lines. 



4.3.2 The Mass Parameter 



In the previous chapter we determined the mass parameters by matching the propagators 
for the zero-frequency modes in the two theories. An equivalent way of determining the 
mass parameter is by demanding that the screening mass in the two theories match. 
Generally, the mass parameter (which is the unphysical screening mass obtained in strict 
perturbation theory) differs from the screening mass obtained in resummed perturbation 
theory, which correctly incorporates the effects of electrostatic screening. We shall need 
the mass parameter m 2 E (A) at next-to-leading order in e. 

In the previous section we saw that the screening mass of the particles is defined as 
the location of the pole of the propagator for spacelike momentum: 

k 2 + n 00 (0, k) = 0, k 2 = -m 2 . (4.17) 

The requirement above implies that 

k 2 + m 2 (A) + U E (k, A) = 0, k 2 = -m 2 , (4.18) 

where He (k, A) is the self-energy of A in the effective theory. One can expand U(k 2 ) = 
n o(0, k) in a Taylor series around k 2 = 0. To determine the screening mass squared to 
order e 4 , we must calculate il'(0) to one loop order and 11(0) to two loop order, and the 
screening mass squared is then given by pSl 



m2«ni(o) + n 2 (o)-ni(o)ni(o), 



(4.19) 
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where Tl n (k) denotes the n'th order contribution to H(k) in the loop expansion. The 
symbol ~ indicates that Eq. (|4.19| ) only holds in strict perturbation theory. The self- 
energy to one-loop order in the full theory reads 



6 ^P} Tr [ P 2 (P + K) 2 
-4(d - 2)e 2 t -L + |(d - 2)e 2 fc 2 i ij + 0(fc 4 /T 2 ). (4.20) 



The corresponding Feynman diagram is shown in Fig. |4.2j . The sum-integrals in Eq. ([4.20 ) 
are standard and are listed in appendix A. The second sum-integral is ultraviolet divergent 
and this divergence may be sidestepped by renormalizing the wave function in the usual 
way. The field strength renormalization constant to the order required is 



This yields 



ri^ 2 ) 
ni(o) 



i 



e 2 T 2 e 2 fc 2 , A 
— — + -^(21n 



12ttV 



12tt 2 
A 



12tt 2 



2 In- 



4ttT 



4ttT 

+ 2 lE - 1 +41n2). 



2 7i; -i + 41n2) + e>(£;7T 2 



(4.21) 



(4.22) 
(4.23) 



Figure 4.2: One-loop self-energy graph in full QED. 

The two-loop expression for the self energy at zero external momentum can be found 
either by a direct computation of the two-loop graphs (See Fig |4.3j ) or by applying the 
formula (see Ref. [ [45|j ) 

11(0) = -e 2 0, (4.24) 

where P is the pressure and \x is the chemical potential. This requires the calculation of 
the free energy to two-loop order including the chemical potential. This has been done 



using contour integration in Ref. [S3], and the two-loop part is 



Po 
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'rrrii i 1° 2rp2 , ^ 
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7T- 



Ti' 



(4.25) 
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We then find 



n 2 (o) 



e 4 T 2 



(4.26) 



Let us just for the record also list the two-loop graphs which are shown in Fig. fO. 
It is interesting to note that diagrams are individually quite complicated, but the sum is 
surprisingly simple: 



II 2 (0) = 4(d-2)e 



{PQ} 



1 i .,. „. 4 ^ 



p 2 q 6 p 2 q 4 



4(rf-2)e 4 J 
Jf 



P{Q} 



1 



p 2 q 6 p 2 q 4 



. (4.27) 



The result is finite before renormalization, and the reason is that the two counterterm 
diagrams cancel as a consequence of the Ward identity (See Fig. fOj). 






Figure 4.3: Two-loop self-energy graphs for n o at zero external momentum. 



Figure 4.4: One-loop counterterm diagrams which cancel due to the Ward identity. 



Now, strict perturbation theory in the effective theory means that the mass term 
should be treated as a perturbation. The corresponding contributions to He{p, A) are 
shown in Fig. 4J3, where the blob indicates a mass insertion. The one- loop diagram 
vanishes in dimensional regularization for massless fields, since the external momentum 
provides the only mass scale in the integral. The matching relation then simply becomes 



mi (A) 



mz 



The mass parameter squared to order e 4 reads: 



m|(A) = T 2 



e 
3~ 



367T 2 



(2 In 



A 
4^T 



+ 2 7E -l + 41n2)- 



8tt 2 



(4.28) 



At this point some comments are in order. Firstly, one could obtain the mass parame- 
ter without renormalizing the wave function as an intermediate step. Instead one uses 
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Eq. (|4.19|) directly and the divergence there is then cancelled by the charge renormaliza- 
tion counterterm. Secondly, by using the renormalization group equation for the coupling 
constant, 

de 2 e 4 

>% = (429) 

one can easily demonstrate that Eq. Q4.28 ) is independent of A. Thus, up to corrections 
of order e 6 , we can replace A by an arbitrary renormalization scale fi. It is interesting to 
note that the mass parameter is equal to the physical screening mass at order e 4 . This is 



in contrast with SQED (see subsection [4.7.2 



Figure 4.5: One-loop self-energy correction in the effective theory. 



4.3.3 The Coefficient of the Unit Operator 



We shall now compute /e(A) to order e 4 in strict perturbation theory and we shall do so 
by matching calculations of In Z in the full theory and in the effective theory. From Eqs. 
|) and ( f4.5|) , we see that the matching condition reads 



hiZ= -f E (A)V + lnZ EQED . 



(4.30) 



Here, Z E q EE) is the partition function of the effective theory, which is given by the path 
integral in Eq. (|4.5|) The calculation of In Z in the full theory involves one-loop, two-loop 
and three-loop diagrams, and we shall discuss them separately in the following. The 



one-loop contribution is displayed in Fig. 4.6, and reads 



--(d-2)t \nP 2 + 2t \nP 2 

2 Tp T{p} 



ll7T 2 T 4 

180 



(4.31) 



Note that the contribution from the ghost field cancels the contribution from two of the 
four polarization states of the photon, and we are left with the contribution from the two 
transverse (or physical) polarization states. The one-loop result is the standard one for a 
gas of noninteracting photons and fermions at temperature T. 



At the two- loop level there is only one vacuum graph (Fig. |4.7[) , which yields 

f . jv — ^/(n = (<i -Dew: - 1 u 1 

r {P} 



tin Tl lT^ = {i ~ 2 ^H lQ} ~ tm ^ 



(4.32) 
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Figure 4.6: One-loop vacuum diagrams in QED. 



Here, £/(-P) is the fermion self-energy function: 



f{Q} Q 2 (P + Q) 2 ' 



(4.33) 



It is interesting to note that this is (d — 2) times the contribution of the corresponding 
diagram in Yukawa theory. The reason is that the fermion self-energy function in QED 
is (d — 2) times the fermion self-energy function in Yukawa theory. 




Figure 4.7: Two- loop vacuum diagrams in QED. 



The three-loop diagrams are displayed in Fig [4.8| , and there are several comments 
that we wish to make; The sum of the first two diagrams in Fig. [4.8| is ultraviolet finite. 



The third diagram in Fig. [|]8] has a linear infrared divergence, which is set to zero in 
dimensional regularization. This diagram is the first in the infinite series of infrared 
divergent diagrams (ring diagrams) that are summed to give the first non-analytic (e 3 ) 
contribution to the free energy. 

The first diagram can also be written in terms of the fermion self-energy function, 
exactly as in the Yukawa case. Recalling the relation between the two fermion self-energy 
functions, the Yukawa result translates into: 



2T {P} 



Tr 



p2 



jpq{k} P 2 Q 2 K 2 (P + Q + KY 
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+2(rf-2)Vi 



QK 



T{p}qk P 2 Q 2 K 2 (P + Q) 2 (P + K) 2 
The second diagram cannot be written in any simple way 



.(4.34) 



4 %{qk} LP 2 Q 2 K 2 (P-Q) 2 (Q-K) 2 (P-K) 2 

( rf -2)(6-rf)e 4 ^ 1 

2 T{pqk} P 2 Q 2 K 2 (P + Q + K) 2 

^ 1 



+(d-2)(d-4)e 



(4.35) 



5pq{x } P 2 Q 2 K 2 (P + Q + K) 2 ' 
The third diagram can also be written in a compact way, involving the polarization tensor: 



ATp P 4 



1 r n2 

n^(P) 



4(d-4)e 



P{Q^} P A Q 2 K 2 

*4 



+i - d ^ e t { PQK}P 2 Q 2 K 2 (P + Q + K) 2 
Tp{qk} P^Q 2 K 2 {P + Q) 2 {P + K) 2 ' 



+16e / 



(4.36) 



Now, let us consider the two-loop counterterm diagrams. These appear in Fig. |4T9| . Gauge 
invariance implies the Ward identity which ensures that Z\ = Z 2 . The correspond- 
ing graphs then cancel, and we are left with the diagram with a photon wave function 
counterterm insertion. This diagram equals the two-loop diagram times (1 — Za). The 
divergence here cancels divergence from the three-loop diagram. Alternatively, one can 
carry out charge renormalization by the substitution e 2 — > Z e 2e 2 in the two- loop graph, 
where 



Z e 2 



1 + 



127r 2 e 



(4.37) 






Figure 4.8: Three-loop vacuum diagrams contributing to three free energy in QED. 
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Collecting our results, we find 
TlnZ 11tt 2 T 4 5e 2 T 4 



V 



T\2 



180 288 16 
8C'(-3) 16 C'(— 1 



-(±- 

7T 2 V 12 



3 V 4vrT' 



- 



319 208 , 

+ In 2 

12 5 



3C(-3) 3 C(-l) 



(4.38) 



The scale A may be traded for an arbitrary scale fi by using the renormalization group 
equation for the running gauge coupling. 






Figure 4.9: Two-loop counterterm diagrams in QED. The first two cancel since Z\ = Z 2 . 

We now turn to the effective theory. The mass parameter is viewed as a perturba- 
tion in the effective theory, as explained above. This implies that lnZ E Q ED is given by 
ordinary one and two-loop diagrams as well as one-loop diagrams with mass insertions 
(which is indicated by a blob in Fig. 4.10|) . The computation is rather simple since loop 
diagrams involving massless fields vanish identically in dimensional regularization. There- 
fore, In 2eqed vanishes in strict perturbation theory and the matching condition turns out 
to be 

rji i 

*-f E {A)T. (4.39) 



V 



/e(A) is then given by minus the right hand side of Eq. ( [4. 38 ) divided by T. 



With the comments after Eq. ([4.38 ) in mind, it is clear that /b(A) has no dependence 
on A at the order we are calculating. The function F = /g(A)T can be viewed as the 
contribution to the free energy from the momentum scale T, which is a typical momentum 
of a particle in the plasma. 
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Figure 4.10: Loop diagrams in the effective theory. 

4.4 The Free Energy and the Electric Screening Mass 

Now that we have determined the short distance coefficients A^(A), m 2 E (A) and /#(A) in 
the effective theory, we can calculate the screening mass squared as well as the free energy 
in QED to order e 5 . In order to do so, we must take properly into account the effects of 
electric screening. This corresponds to the following decomposition of the Lagrangian: 

(£eqed)o = -F ij F ij + -(d i A y + -m 2 E (A)A 2 + £ g t+£ gh (4.40) 

(£EQED)int = ^-At + SC. (4.41) 

The electric screening mass is given by the location of the pole in the propagator and 
at leading order it is simply m s = m#(A). To leading order, the self-energy function 



He (A;, A) is given by the first Feynman diagram in figure [4.5| and is independent of the 
external momentum. Eq. ( |4.18| ) then turns into 

m\ = m 2 E + U E (k,A) 

= ml + ^[ — (4.42) 
2 J P p z + m E 

Using this result, the expression of Ag(A), and expanding the mass parameter in powers 
of e, we obtain the electric screening mass squared to order e 5 : 

< = T " [t ~ 4~2 ( 2 ln Tt + 2lE ~ 1 + 4 ln 2) ~ T* + 7^1 • ^ 443 ) 
L 3 36^^ 4nT 87H 4\/37r 3J 
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It is easily checked that the result is RG-invariant, as required. Furthermore, our result 
agrees with the calculation of Blaizot et al ]52]. Note also that there is no e 3 term in the 
expression for the screening mass squared in contrast with both 4 -theory and SQED. 
The reason is that there are no bosonic propagators in the one-loop self-energy graph 
in QED and fermions need no resummation, since their Matsubara frequencies are never 
zero. Finally, the e 5 term is non-analytic in e 2 and so corresponds to the summation 
of an infinite number of diagrams in terms of bare perturbation theory. The diagrams 
are the two-loop graphs in Fig. |4.3j with insertions of any number of il o(fco;^) on the 
internal photon lines. This is in complete analogy with the infinite string of diagrams 
which we discussed in connection with the evaluation of the screening mass in Yukawa 
theory in chapter two. This is a good example of the efficiency of the effective field 
theory approach. Instead of summing an infinite number of diagrams, we simply perform 
a one-loop computation in three dimensions. 



The calculation of the free energy in the effective theory is straightforward. Bearing 
in mind the fact that the self interaction term contributes to the free energy first at order 
e 6 , we only need to perform a one- loop computation, and so we obtain 



T In Z E qed 
V 



m, 



■(d-3)T f hip 2 . (4.44) 



The relevant Feynman graphs are displayed in Fig. |4.10|, except that the diagrams with 
a mass insertion, as well as the two-loop graph are not included. The contributions from 
the gauge field and ghost vanish. Using the expression for the mass of the scalar field and 
expanding it in powers of e yields the following contribution to the free energy: 



TlnZ E QED _ e 3 T 4 



V 



(4 In 



A 
4^T 



+ 4 7i; + 7 + 81n2 



(4.45) 



This term takes into account the effects from long distance scales of order l/(eT), which 
can be associated with the scale of electric screening. Using Eqs. ( |4.30|) , (|4.83| ) and ( |4.45| ), 
one finally obtains 



T\nZ 
V 



11tt 2 T 4 5e 2 T 4 e 3 T 4 



180 
319 
' ~12 



288 



+ 



36^ 
8C'(-3) 



In 2 + 

5 3 C(-3) 



16tt 2 V 12 / 

i6 C(-i: 

3 C(-i) 



3 MttT' 



e 5 T 4 



576^3 



4 In 



- 47^ 
A 



4ttT 



+ 4 7 £ + 7 + 81n2 



(4.46) 



This result is renormalization group invariant as required for a physical quantity. This can 
easily be checked by using the one-loop /3-function in QED. Moreover, it is in agreement 
with the computation of Parwani f92|| , and Zhai and Kastening who use resummed 
perturbation theory. The advantage of the effective field theory approach should now 
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be clear; In order to to extract the e 3 and e 5 contributions to the free energy using 
resummation, one must use the resummed propagator in every diagram and find the 
subleading pieces by subtracting the leading ones, as we indicated in the previous chapter. 
This is at least a rather tedious task. Here, we obtain the non-analytic terms in the free 
energy from a straightforward one-loop calculation. 

It is interesting to note that there are no terms in the expressions for the screen- 
ing mass or the free energy which involve logarithms of the coupling constant. This is 
in contrast with QCD and SQED, where the free energy contains a term proportional 
to g 4 lng [f75fl . This can be understood in terms of the renormalization group and the 



renormalization of the parameters in the effective theory. 

The short-distance coefficients in £eqed are obtained by integrating out the non-zero 
Matsubara frequencies, and they are polynomials in ln(A/47rT). So in order to avoid large 
logs, we must choose the cutoff of order T or 2ixT . The latter is a more physical choice 
in the sense that the heavy modes have masses 2ttT or more. 

Moreover, in the effective theory one generally encounters logs of A/mg in the per- 
turbative expansion. This implies that we must choose the scale A in the effective theory 
of order to control the perturbative expansion. Thus, one must take the parameters 
in EQED from the scale 1nT to the scale the, using the equations which govern their 
evolution with the scale. These evolution equations follow from the requirement that 
physical quantities be independent of the cutoff, and are in the form 

= &(C(A)). (4.47) 

The beta-functions can be written as power series in the coupling constants of the effective 
theory. Using dimensional arguments, we can infer the general structure. Consider first 
m E (A). Its expansion must be a quadratic polynomial in e E (A) and A#(A), and other 
coefficients, so that the dimension of every term is two. The only coefficient that contribute 
to order e 4 is ejj(A) (Recall that As (A) ~ e 4 and that other coefficients are of even higher 
powers of e). Reasoning along the same lines, we conclude that the beta-function of /e(A) 
involves the terms e|*(A)m|;(A), A^(A)m|;(A), as well as cubic polynomials in e%(A) and 
Xe{A) (and other terms which have dimension three). The only relevant term at the 
order e 4 is e|;(A)m|;(A). Finally, we mention that the beta-functions for eg(A) and Xe{A) 
vanishes for superrenormalizable interactions |75| . This implies that the beta- functions 



are highly suppressed by powers of the coupling, since these receive contributions only 
from higher order operators. 

However, we have already noted that m E (A), /e(A) and A^(A) are independent of A 
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at this order so we have 

A*^ = £>(«•), A^=0(e«), A^ = 0(A (4.48) 

The dependence of A in the parameters in £eqed are canceled by the cutoff dependence 
in the effective theory, and so the vanishing of /3-functions for these parameters explains 
that no logs of e occur. More specifically, the fact that m E (A) does not run implies 
the non-existence of a term e 4 lne. Moreover, since m 2 E (A) has an expansion in e 2 and 
(A) does not run, there can be no term e 5 lne term either. Similarly, the vanishing of 
/3-function for the unit operator is responsible for the fact that a term e 4 In e is absent in 
the expression for the free energy. Since /75(A) has an expansion in even powers of e, an 
e 5 lne can only arise from the mass parameter. However, since m E (A) does not run at 
next-to-leading order in e, this explain the absence of such a term. 



4.5 QED Versus QCD 

Finally, we would like to discuss a computational as well a principal difference between 
QED and QCD, when calculating the free energy beyond the fifth order in the coupling 
constant. In QCD the computation of the free energy involves the construction of two 
effective field theories, which reflects the fact that there are contributions from three 
different momentum scales (T, gT and g 2 T, where g is the gauge coupling) [75]. The first 
effective field theory, called electrostatic QCD (EQCD), consists of the magnetostatic 
Al field and the electrostatic field Aq. The unit operator /73(A) as well as the other 
parameters in EQCD are then determined by the usual matching procedure and /75(A) 
gives the contribution to the free energy from the short distance scale 1/T. 

The second effective field theory is called magnetostatic QCD (MQCD) and consists 
simply of the self-interacting magnetostatic gauge field A°:. Again, the unit operator 
/m(A) and the coupling constants of this effective theory can be determined by matching 
calculations, and /m(A) yields the contribution to the free energy from the distance scale 
l/(gT). Now, the perturbative expansion in MQCD is plagued with infrared divergences, 
implying that the functional integral can only be calculated non-perturbatively, e.g. by 
putting MQCD on a lattice. Using lattice simulations the path integral may be computed, 
so that one obtains the contribution to the free energy from the scale l/(g 2 T). It can be 
written as a power series in g starting at order g 6 . The leading contribution, 0(g 6 ), has 



very recently been determined by Karsch et al. p6 \. 



One can, of course, construct a second effective field theory, which naturally is termed 
magnetostatic QED (MQED), but it is completely unnecessary. Although this is obvious 
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from a physical point of view (there are only two scales in QED), it is instructive to see 
this in practice. Now, MQED contains all operators which can be constructed out of the 
fields Ai, which satisfy the symmetries, such as gauge invariance and rotational symmetry 
. We can then schematically write 



Aviqed — + 9M(-^)FijV 2 Fij + h M (A) (FijFij) 2 + ... 



(4.49) 



The second operator is the analog to the Uehling term which is well-known from non- 
relativistic atomic physics |97[|. The coefficient has be determined by Landsman |7I] and 
is 7C(3)e 2 /9607r 4 T 2 . The third term corresponds to one of the operators in the famous 
Euler-Heisenberg Lagrangian P3|. Its coefficient has been worked out in Ref. P5[ , but it is 
in an extremely complicated form. Nevertheless, we shall consider it as an example. The 



two-loop graph arising from this interaction is shown in Fig. [4.11] . Using power counting 
arguments it is easy to verify that only power ultraviolet divergences occur. The canonical 
dimensionality of is 3/2, and so h M (A) must be proportional to e 4 /T 3 (at leading order 
in e). The contribution then goes like 



fip, q) 

T 3 J pq p 2 q 2 



(4.50) 



Here, f{p,q) has dimension four, since the dimension for the free energy is three. The 
power ultraviolet divergences are artifacts of the regulator, and in dimensional regular- 
ization they are set to zero p5[ . Thus, the contribution vanishes. One can use simi- 
lar arguments to conclude that every operator, except for the unit operator, gives zero 
contribution to the free energy in magnetostatic QED. The point is now that we are 
actually computing the unit operator, /m(A) when we do perturbative calculations in 
EQED and include the mass term m£(A) in the unperturbed part of the Lagrangian. 
Thus, there is no need for determining the other coefficients in MQED. We then have 
F = f E (A)T + f M (A)T. 




Figure 4.11: Example of a vanishing two-loop diagrams in magnetostatic QED. 
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Let us close this section by some comments on the infrared catastrophe in QCD and 



its solution by the present method [75]. It is a well-known fact that the free energy of 
nonabelian gauge theories may be calculated to fifth order in the coupling using resummed 
perturbation theory. However, the method breaks down at order g 6 , due to infrared 
divergences, as first pointed out by Linde p9[ . These divergences arise from regions 
where all internal energies vanish, and so the singularities are the same as in three- 
dimensional pure QCD. Thus, the breakdown of perturbation theory simply reflects the 
infrared problems appearing in a perturbative treatment of any nonabelian gauge theory 
in three dimensions, in particular MQCD (although it is well behaved nonperturbatively 
with a mass gap of order g 2 T). In the present approach one can compute order by order the 
contributions to the free energy, although some coefficients must be evaluated numerically. 
The infrared problems can naturally be avoided if one uses lattice simulations directly in 
four dimensions. However, this is extremely time consuming in comparison with MQCD, 
and the time savings here arise from the reduction of the problem from four to three 
dimensions, and also by integrating out the fermions. 



4.6 Electrostatic Scalar Electrodynamics 

In this section we continue our study of effective field theories by investigating scalar 
electrodynamics. The Euclidean Lagrangian of massless SQED is 

£sqed = \f^F, u + (P^)t(P^) + ^($t$) 2 + £gf + £ gh , (4.51) 

The effective field theory is called electrostatic scalar electrodynamics (ESQED), and 
consists of a real massive scalar field (the temporal component of the gauge field) coupled 
to scalar electrodynamics in three dimensions. According to the preceding discussion we 
must write down the most general Lagrangian which respects the symmetries at high T. 
Lorentz invariance is broken at finite temperature, so we must allow for a mass term 
and self-interactions for the timelike component of the gauge field. Moreover, the Ward 
identity at high temperature implies that the effective Lagrangian is a gauge invariant 
function of the fields A{. Finally, there is a ^-symmetry for the fields and Aq. This 
symmetry follows from the corresponding symmetries in the full theory. The effective 
Lagrangian then has the general form 

£esqed = \F i3 F l3 + {V4)\V4) + M\X)^<P + ^{d l A a ){d l A a ) + 

l -m 2 E {A)Al + e|(A)0V4j + ^(0V) 2 + C gi + C gh + SC. (4.52) 
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Here SC represents all other terms consistent with the symmetries. Examples of such 
terms are \e(A.)Aq, which is superrenormalizable and /i^(A)(FjjFjj) 2 , which is nonrenor- 
malizable. 



4.7 The Parameters in ESQED 

The coefficients in ESQED are again determined by ordinary perturbation theory in pow- 
ers of e 2 and A, neglecting resummation. In the full theory we then split the Lagrangian 
into a free part and an interaction part accordingly: 

(£sqed)o = + (^$)t(^$)+£ gf + £ gh , 

(£s Q ED)int = e2$ t$ K - ieA^&dfi - <S>d^) + ^) 2 . 

In the effective theory the masses as well as higher order operators are treated as pertur- 
bations. We then write £esqed = (£esqed)o + (^ESQED)int an d strict perturbation theory 
corresponds to the following partition of the effective Lagrangian 

(£esqed)o = \ F a F ij + W) f W) + ^(d t A ) 2 + £ gf + £ gh , 
(^ESQED)int = \rn 2 E (A)J% + M 2 (A)0V + e|(A)0V(A? + Al) 

-ie E (A)M^9i(l> - <t>ditf) + ^V<£) 2 + SC. (4.53) 

o 

In SQED the scalar field is denoted by a dashed line, the photon by a wiggly line, and 
the ghost by a dotted line. In ESQED we have the additional convention that the real 
scalar field is indicated by a solid line. 



4.7.1 The Coupling Constants 



For the present calculations, we need the gauge coupling eg (A) only to leading order in e 
and A. By using the relation between the gauge fields in the two theories 




(4.54) 



and comparing £esqed with JJf dr Csqed, we find 



4(A) = e 2 T. 



(4.55) 
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At this order there is no dependence on the renormalization scale A. Similarly one finds 
at leading order 

A 3 (A) = AT. (4.56) 



4.7.2 The Mass Parameters 

In this subsection we calculate the parameters M 2 (A) and m 2 E (A) at leading and next-to- 
leading order in the coupling constants e and A, respectively. The physical interpretation 
of a mass parameter is that it is the contribution to the physical screening mass from 
momenta of order T. The simplest way of determining the mass parameters is to match 
the screening masses in SQED and in ESQED. Denoting the self-energy for the field $ by 
E(/c ,k), the scalar screening mass is the solution to the equation^ 

£; 2 + E(0,k) = 0, k 2 = -m 2 . (4.57) 

The matching requirement implies that 

£; 2 + M 2 (A) + E(A;,A) =0, k 2 = -m 2 , (4.58) 

where £ (k, A) is the self-energy of the field in the effective theory. The self-consistent 
solution to Eq. (|4.57|) is to leading order in the coupling constants m 2 ~ £i(0). Here 
E n (fc 2 ) = £ n (0, k) denotes the nth order contribution to the self-energy in the loop ex- 
pansion and the symbol ~ is a reminder that this unphysical screening mass is obtained 



in strict perturbation theory. The relevant diagrams are depicted in Fig. |4.12| and the 
one-loop self-energy at zero external momentum is given by 

E 1 (0) = (rf-l)e 2 ^^ + ^^. (4.59) 

The limit d — > 4 is perfectly finite, and this immediately gives 

= , x e 2 T 2 AT 2 , , 

Si(0) = ^ + ^. (4.60) 

The self-energy function A) vanishes in strict perturbation theory, since all the 
propagators are massless. Hence the matching requirement gives m 2 « M 2 (A) and the 
mass parameter is 

e 2 T 2 XT 2 

M 2 (A) = e -L + ±L-. (4.61) 



2 We remind the reader that this screening mass has nothing to do with the screening of electric fields. 
As previously noted, this is a quantity which gives information about the screening of static scalar fields 
due to rearrangements in the plasma. 
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Figure 4.12: One-loop scalar self-energy diagrams in SQED. 

At this order M 2 (A) is independent of the scale A. 

Let us now turn to the mass parameter m 2 E (A). The screening mass is again defined 
as the pole of the propagator at spacelike momentum 



A; 2 + n 00 (0,k) = 0, 



-m 2 s . 



(4.62) 



The self-energy function is given by a series expansion in e 2 and can also be expanded in a 
Taylor series around k 2 = 0. The self-consistent solution to Eq. (|4.62|) at next-to-leading 
order in the coupling constant is in analogy with QED 



m 2 ~ 



n x (o) + n 2 (o) i-ni(o) 



(4.63) 



Here, we have again defined U(k 2 ) = IIoo(0, k) and U n (k 2 ) denotes the nth order contri- 
bution to I1(A; 2 ) in the loop expansion. The one-loop self-energy is shown in Fig |4.13| . and 
equals 



Mk 2 



2e' 



P 2 



-ke- 



p p 2 (p + ky 



(4.64) 



Expanding in powers of the external momentum and integrating by parts in d — 1 dimen- 
sions yields 

Mk 2 ) = ^%j- 2 -^%§i + ^ 2 k%§- 6 +0(kVT 2 ). (4.65) 

The last sum-integral is ultraviolet divergent and this divergence may be removed by the 
wave function renormalization counterterm: 



Z A = l~ 



3(4tt)V 



One then obtains 

n;(o) 



e 2 T 2 



+ 



2A; 2 



3(4tt) : 



;(ln 



A 
4^T 



+ lE + l) + 0(k 4 /T 2 



3(4vr 



2 /, A 

(In— — + JE + 1)- 



AttT 



(4.66) 



(4.67) 
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Figure 4.13: One-loop self-energy diagrams in the full theory. 

We also need the self-energy at zero external momentum to two loop order. The 
contributing diagrams are displayed in Fig. |4.14j . Many of the two-loop sum-integrals 



vanish in dimensional regularization, while others factorize into products of one-loop sum- 
integrals. After some calculations we find 

n 2 (0) = 8(rf-l)e 4 i _gL-2(d-l)e 4 i 



pqP 6 Q 2 y ' TpqP 4 Q 2 

4Ae 2 r 1 16Ae 2 r pi , . 

+ ^^t 75§^- ( 4 - 68 ) 



3 TpqP^Q 2 3 TpqP^Q 2 

The ultraviolet divergences in the above sum-integrals actually cancel, and so we are 
left with a finite expression for 112(0). This cancelation is exactly the same as the one 
we encountered in QED and it reflects the Ward identity. Using the tabulated one-loop 
sum-integrals in appendix A, one obtains 

e 4 T 2 Ae 2 T 2 

n 2 (0) = ^~ 2 + (4.69) 
(4tt) z 727T 

Using these results, we finally obtain m 2 to order e 4 : 

e 2 T 2 r /2 , A 2 7x e 2 t Ae 2 T 2 



m 2 



1 -(3 ln 4^T + 3^-3)(4^] + ^- ( 4 ' 7 °) 



72tt 2 



Note that one could have obtained this result without carrying out wave function renor- 
malization in Eq. ( |4.64| ). Instead one uses Eq. ( |4.63| ) and the divergence there is canceled 
by the charge renormalization counterterm. 

In the effective theory the contributing diagrams are the usual one - and two-loop 
graphs plus one-loop graphs with mass insertions. Denoting the self-energy by n#(A;, A), 
the screening mass is given by the solution to the equation 

k 2 + m|(A) + U E (k, A) = 0, k 2 = -m 2 . (4.71) 



We are now familiar with the fact that all loop integrals involve massless fields and these 
vanish in dimensional regularization. Hence Yi^{k,K) = 0, and so the matching relation 
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becomes m 2 (A) « m 2 . Thus 



m|(A) 



e 2 T 2 r 



2 _A_ 2 7x e 2 - 



2rr2 



+ 



Xe T 
72tt 2 



(4.72) 



One can verify that the apparent A- dependence of m 2 E (A) is illusory. This implies that, 
up to correction of order e 6 , one can trade A for an arbitrary renormalization scale \i. The 
reason behind this fact is that the physical screening mass does not receive logarithmic 
corrections in the effective theory to order e 4 (see section [O) . 





Figure 4.14: Two- loop self-energy diagrams in the full theory. 
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4.7.3 The Coefficient of the Unit Operator 



In this subsection we shall determine the coefficient of the unit operator. We shall consider 
the one, two and three-loop contributions as well the contributions from the counterterms 
diagrams separately. The matching condition we use to determine /e(A) follows from 
the two path integral representations of the partition function in complete analogy with 
QED: 

InZ = -f E (A)V + lnZ ESQ ED. (4.73) 
Let us first focus on the SQED. The one-loop graphs in the underlying theory are depicted 



in Fig. |4.15| and the corresponding contribution reads 

2tt 2 



45 



(4.74) 




Figure 4.15: One-loop vacuum diagrams in SQED. 
The two-loop diagrams are shown Fig. |4.16| . After some purely algebraic manipula- 



tions they factorize into products of simpler one-loop sum-integrals. The result is 

,2. 



ztp Q P 2 Q 2 + ^ 2 )e2 ^ Q P 2 Q 2 



/TV-' -A 5r~ 
U2/ 



3 + 2 



(4.75) 



Let us now turn to the three-loop diagrams. These are displayed in Fig. |4.17| . Here, 
the shaded blob means insertion of the one-loop polarization tensor U fll/ (ko, k), while the 
black blob implies insertion of the scalar self-energy function S(fc ,k), also at one loop. 

The first four diagrams can expressed entirely in terms of the bosonic basketball. 
After some purely algebraic manipulations one finds 



+ (d - 13/4)e 



The fifth diagram gives a contribution 



ATp P 4 



1 r n2 
IWP) 



pqk P 2 Q 2 K 2 (P + Q + K) 2 



Tpqk P 4 Q 2 K 2 (P - K) 2 (P - Q) 2 



(4.76) 



1 



4 Tpqk P 2 Q 2 K 2 (P + Q + K) 2 ' 



pqk P A Q 2 K 2 
(4.77) 
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Figure 4.16: Two-loop diagrams for the free energy. 



The last graph reads 



Tpqk P 2 Q 2 K 2 (P + Q + K) 2 



2A 2 2(d-l)\e 2 (rf-l)V- 



9 



+ 



+ 



pqk P 4 Q 2 K 2 



.(4.78) 



/ 

/ 

\ 

\ 



\ 
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/ 

/ 




1 t^^^O 
1 ^^^^7 




• 



Figure 4.17: Three- loop diagrams for the free energy. 

With our experience of both Yukawa theory and QED in mind, we carry out renor- 
malization by the substitutions A — > Z x \ and e 2 — > Z e 2e 2 . To the order needed, the 
renormalization constants are: 



Z\\ = A + 



5A 2 - 18Ae 2 + 54e 4 
487r 2 e : 



1 + 



487rV 



(4.79) 
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After we have carried out coupling constant renormalization, we are still left with a term 
proportional to 1/e. This term is canceled by 6Je{A), which is the counterterm for /#(A). 
This is the only nonvanishing term in the calculation of In ^esqed m strict perturbation 
theory (we remind the reader that this involves massless fields, and the loop-integrals 
are therefore zero). According to Ref. |p8| , <5/s(A) can be computed by considering the 
ultraviolet logarithmic divergences in the effective theory, when ones uses dimensional 
regularization. Generally, 5/^(A), is a power series in M 2 (A), m E (A), e E (A) and A3 (A). 



At leading order it turns out that it is given by 

Sf E (A) - 



e|M 2 1 



2(4vr) 



(4.80) 



which follows from a two-loop calculation in the next section. Since the mass M is 
multiplied by 1/e it is necessary to expand it to first order in e when expressing 5f E {A) 
in terms of e 2 , A and T. From Eq. ( f4.59| ) one finds 



<9M 2 



<9e 



e2T2 f 6 .„. A 



e=0 

This implies that 



12 



4ttT 



+ 4 + 6 



c(-i; 



c(-i; 



AT 2 , 
+ 2 In • 



A 
4+T 



2 + 2 



C'(-i) 



Sf E (A)T 



(4tt) 2 V 12^ 
Ae 2 /T 2 xr4 



^ + 361n^ 

e 4ttT 



+ 24 + 36 



C(-i 
C'(-i 



(4.81) 



C(-i) 



V19 



(4tt) 2 V12^ 

Putting our results together, one finally obtains 



, A 
- + 8 In — - 

e 4vrT 



+ 8 + 



,c'(-i: 



c(-i; 



(4.82) 



f E (A)T 



2vr 2 T 4 



T 2 \ 2 r A 



45 
A 2 



3 + 



5e 2 



T 2 N2 r 10. A 



16?r 2 V 12 ) 
16vr 2 U2^ 



31 2C(-3) 4C'(-i; 
9 4vrT 9 lE 45 3 C(-3) 3 C(-l) 



2r , A 

12 In 



— (—) 

16vr 2 U2/ 



T 2 \ 2 r257 



20 

+ A'Je H 

4ttT 1 3 

A 



>C'( 



In — — + 137r 

3 4ttT ' 



C(-i) 

164 110C'(-3) 



328 C'( 



. (4.83) 



3 3 C(-3) 3 C(~l] 
In contrast to the corresponding three-loop calculation in QED, /e(A) is not renormal- 
ization group invariant. This follows easily from the renormalization group equations, 
and the reason is that there is a logarithmic ultraviolet divergence at two-loop order in 
ESQED. 



de 2 
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Instead, /e(A) satisfies an evolution equation ffSH , which follows from Eq. ( |4.80|) 



d f E (A) _e|M 2 
dk 2(4vr) 



^ ja = -^i- (4-86) 



4.8 Calculations in ESQED 

Now that we have determined the short-distance coefficients we shall use the effective 
three-dimensional field theory and calculate the electric screening mass and the free energy. 
We shall do so using perturbation theory and in order to take the physical effect of 
screening into account, we must again include the mass parameters in the free part of the 
effective Lagrangian. This corresponds to the following partition of £esqed : 

(£esqed)o = ^F iJ + (9 i 0t)(^0) + M 2 (A)0t0 + I(^ o ) 2 

+ ^m 2 E (A)Al + £ gi + £ gh , (4.87) 

(^ESQED)int = e 2 E (A)^4>(A 2 + A 2 ) - le E {k) A^ - 4>d^) 

+^S-^^ + 5C. (4.88) 
o 



The physical screening masses are given by the self-consistent solutions to Eqs. ( f4.58| ) 



and ( |4.71| ). The solution to Eq. ( f4.58| ) to leading order in coupling is equal to the mass 
parameter M 2 (A). However, recently it has been realized that this equation has no self- 
consistent solution beyond leading order in perturbation theory [52,53]. The problem is 
the last diagram in Fig. |4.12| , which has a branch point singularity at k — im s . The 
problem is the same as in QCD, namely a scalar field interacting with a massless gauge 
field in three dimensions. In QCD this singularity may be screened by a magnetic mass 
of nonperturbative origin. In SQED the magnetic mass is absent since it is an Abelian 
theory f93fl , and so the problem cannot be solved this way. We shall not discuss this any 



further, but refer to Ref. |53| where a nonperturbative definition of the scalar screening 
mass is discussed in detail. 

The one and two-loop diagrams that contribute to the electric screening mass in 



ESQED are displayed in Figs. |4.18| . We then find 



1 _ 4 r Si. 

E 



\2 



p 2 + M 2 Jpq q 2 (p 2 + M 2 ) 



+2ei [ \» ^ _ -2e1 



-4e| 



pq (p - q) 2 {p 2 + M 2 ) 2 {q 2 + M 2 ) * Jpq {q 2 + M 2 ) 2 {p 2 + m% 



pq {p 2 + M 2 )(p 2 + q 2 + m|) [(p + q + k) 2 + M 2 } 
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4A 3 e| r 1 



,, q (p 2 + M 2 )(g 2 + M 2 ) 2 ' (4 ' 89) 

The integrals may be reduced to known ones by algebraic manipulations, which involve 
changes of variables. The integrals needed are tabulated in appendix B. The second inte- 
gral above, which corresponds to the the fifth graph in Fig. [4.18| , vanishes in dimensional 
regularization due to the masslessness of the photon. Moreover, the fifth integral is depen- 
dent on the external momentum. This is the same as the scalar setting sun diagram we 
met in the chapter on resummation, albeit with different masses. We recall that the self- 
consistent solution to Eq. (|4.71| ) is found by evaluating the integral at the point k = im s . 
The calculation of this diagram is carried out in some detail in appendix C. Notice also 
that the logarithmic divergence from this integral is exactly canceled by a corresponding 
term in the second two-loop integral above. Adding the different pieces, we obtain the 
physical screening mass squared to order e 4 and Ae 2 : 

-e 2 Me 2 e 4 / m M, , . m N \ 

-3-^ + (2^(- 1 + 4M + (1+ m )k(1+ M ) ) 

2e" , , A m Xe 2 T 2 

: (l +7E + l„ )]+ (4.90) 



ml = T" 



(12vr) 2 V irj 4nT< 



Settting A to zero, our calculations reproduce the result of Blaizot et al. |52| , who used 
resummation methods. The inclusion of a scalar self-interaction term only produces a 
modification of the scalar mass parameter M 2 (A) proportional to A, and an almost trivial 
term in the expression for the electric screening mass squared, proportional to Ae 2 . 

Using the renormalization group equation for e, we find that the physical screening 
mass is independent of the renormalization scale A up to corrections of order e 5 . We have 
for completeness also checked that the incorrect definition m 2 = n o(0,k — > 0) does not 
satisfy the RG-equation, exactly as in QED, first pointed out by Rebhan |51| . The lesson 
we can learn from this, is that gauge fixing independence is only a necessary, but not a 
sufficient criterion for a quantity to be physical (recall that H^ko, k) is manifestly gauge 
fixing independent in Abelian gauge theories). 

Let us now move on to the calculation of the free energy in ESQED. The one and 
two-loop contributions are depicted in Fig. |4.19] and yield 



r '"^ ESQED = -\Tj\^ + ml)-Tjln{j? + Kfi)-\(d-?,)TjW 
4eiT f <P + q) l„ „ - delT 



pq (p 2 + M 2 )(g 2 + M 2 )(p-q) 2 * J pq (p 2 + M 2 )g 2 

A* r l 



pq (p 2 + M 2 )(g 2 + m 2 E ) 3 J pq {p 2 + M 2 )(g 2 + M 2 ) 



-5f E (A). (4.91) 
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Figure 4.18: One and two-loop self-energy diagrams in the effective theory. 



The two-loop contributions may be reduced to products of one-loop integrals and two-loop 
integrals which are tabulated in appendix B. Using this and Eq. ( |4.83 ) we obtain the free 
energy through order A 2 , Ae 2 and e 4 : 



27r 2 r 4 

45 
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+ 
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167T 2 \12' 

Ae 2 /T 2 x2 r 
16tt 2 ^12^ 
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r A 5e 2 i 
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4ttT 



e 3 T 4 
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M 3 T e 3 MT 3 
6tt 16tt 2 V3 
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12 In 



A 



A 52 
- 16 In h 47 S h 

2M /£ 3 



3 C(-l 

,C'(-i) 



16tt 2 v 12 / 



257 A , A 2 

In — - - 72 In — + 137^ + - 

3 AttT 2M ' 3 



C(-i) 
no C'( 



-3) + 328C'( 



3 C(-3) 3 C(-l) 



Firstly, we note that the two-loop contribution in the gauge sector is exactly the same 
as in spinor QED. Secondly, our result is renormalization group invariant up to order 
A 2 , Ae 2 and e 4 , as it must be. This can be easily checked by using the RG-equations 
for the running coupling constants. Thirdly, we notice the appearance of e 4 ln(A/M) 
and Ae 2 ln(A/M) terms. These are necessary in order to cancel the A-dependence in 
/e(A). The fact that logarithms of the coupling constants occur is then attributed to the 
renormalization of /^(A). A corresponding term g 4 ln(A/m£) was also found by Braaten 
and Nieto in the case of QCD [75]. No terms of order A 2 ln(A/m) arise in A0 4 -theory, 
since /(A) does not run at next-to- leading order in A 5^ . 



Moreover, from chapter three, we know that the scalar mass parameter M 2 (A) is not 
renormalization group invariant, in contrast with the corresponding mass parameter in 
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QCD. From the general arguments given in section 14, this explains the absence of a 
g 5 Ing term in expression for the free energy in QCD, and allows us to predict such terms 
in the free energy in SQED (just as in pure A0 4 -theory). 



4.9 Summary and further Outlook 

In this chapter we carried out detailed calculations in gauge theories at finite temperature 
using effective field theory to unravel the contributions to physical quantities from the 
scales T and eT ~ a/AT. This is the advantage of the effective field approach over the 
more conventional resummation procedure; the latter complicates calculations unneces- 
sary because the sum-integrals involve both T and m. The simplifications of effective field 
theory are more transparent as we go to higher orders in the loop expansion. 

In contrast with QCD, one can, in principle, compute the free energy in Abelian gauge 



theories to any order desired in the coupling constant |75|. So let us briefly outline what 
it takes to push the calculational frontier to the next order. Consider first QED. In order 
to obtain the free energy to sixth order, we must know /^(A) to order e 6 . This requires 
the calculation of four- loop diagrams in the full theory. The mass parameter m 2 E (A) 
contribute only at odd powers in e in a one-loop calculation in EQED. This follows from 
the facts that one can write its coefficient in powers of e 2 , and that the free energy at 
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one-loop goes like m\(N)T. So the next contribution is first at order e 7 . Moreover, 
there is only one operator in addition to /e(A) that contributes at order e 6 , namely the 
quartic self-interaction of A . Since we already know its coefficient, we are left with a 
straightforward two- loop calculation in EQED (the double bubble). Hence, the challenge 
is the four-loop diagrams in full QED. 

In the case of SQED, we need the parameters /#(A), m 2 E (h) and M 2 (A) to order 
A 2 , Ae 2 and e 4 . We already know these coefficients, although the latter parameter was 
determined in chapter three, using the effective potential. There are no new operators 
which contribute at this order, so all that remains is calculating the three-loop diagrams 
in ESQED. Moreover, the evolution equations should play a greater role here than in 
QED, since the beta-function of the scalar mass parameter is non-vanishing, exactly as 
in the somewhat simpler A0 4 -theory |5£|. In particular, they can be used to sum up 
leading logarithms. I have not yet performed the three-loop calculation, but it should be 
a manageable task, and will be the subject of further investigation. 
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Chapter 5 



Conclusions 



It is time to summarize and draw some conclusions from the present work. In the first 
chapter, we studied a rather old problem, namely that of charged particles in external 
magnetic fields, to demonstrate calculations where the results could be expressed in closed 
form. The fermions in two spatial dimensions show nontrivial properties such as de-Haas 
van Alphen oscillations, induced vacuum charge and induced gauge-noninvariant terms in 
the effective action. The bosons are somewhat more trivial, except for the fact that they 
go from a diamagnetic phase to a paramagnetic one. 

The second chapter was devoted to resummation, which has been the dominant way of 
doing consistent calculations at finite temperature. We applied the formalism to calculate 
the screening mass and the free energy in Yukawa theory, which are both static quantities. 
In many respects, resummation has been beaten by effective field theory, but it will be 
with us for years to come. It it still the only way of solving dynamical problems at 
finite temperature, since these rely on an unambiguous analytic continuation of Greens 
functions from imaginary time to real time, or direct calculations in real time. 

The bulk of the thesis is on effective field theory. I have tried to bring about the philos- 
ophy and the general ideas about effective field theory by some examples and discussion. 
The effective field theory program has led to an increased understanding of quantum field 
theory, conceptually speaking, and a better insight into concrete physical systems. This 
approach opens up a new way of attacking difficult problems in physics. I am sure that 
the future will provide examples where effective field theory either solves problems which 
cannot be solved by more traditional methods, or is used to push the calculational frontier 
to higher orders in perturbation theory. The message is (at least) two-fold: 

Firstly, there is the modern view on renormalization which is intimately connected 
with the fact the every quantum field theory should be looked upon as an effective de- 
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scription at some scale. From this perspective renormalizability is no longer a requirement 
of a consistent and useful quantum field theory. This has given meaning to gravity and 
the nonlinear sigma model beyond the tree approximation. The next generation of text 
books certainly should be rewritten on this point! 

Secondly, I have tried to demonstrate the efficiency of the effective field theory pro- 
gram: identify the scales in your problem and take care of them, one by one, by integrating 
them out successively. Thus, effective field theory unravels the contributions to physical 
quantities from different momentum scales, and has become an important tool in prac- 
tical calculations. In order for this program to work, the energy scales in the system 
must be widely separated. I have applied these ideas to quantum field theories at finite 
temperature, and reproduced known results relatively easily. I have also obtained new 
results with the effective field theory approach which minimizes the calculational efforts. 
Effective field theories are here to stay! 



Appendix A 

Sum-integrals in the Full Theory 



In this appendix we summarize our conventions and define the sum-integrals used in 
the calculations. We use the imaginary time formalism, where the four-momentum is 
P = (po, p) with P 2 = pi + p 2 . The Euclidean energy takes on discrete values, po = 2mrT 
for bosons and po = (2n + l)irT for fermions. Dimensional regularization is used to regu- 
larize both infrared and ultraviolet divergences by working in d = 4 — 2e dimensions, and 
we apply the MS renormalization scheme. 

The shorthand notations for the sum-integrals in the bosonic and fermionic cases are, 
respectively: 



po=2irnT 



m = C^)'t e f^nn (a. 2) 



t P > V 47F 1 P0 =2,{n+1/2)T J i 2 *) 

All of the sum-integrals used in this thesis have been calculated and tabulated by Arnold 



and Zhai [42|, Zhai and Kastening [44|, and Braaten and Nieto [75]. 



The general formulas for the bosonic one-loop sum-integrals are 
1 ^VvT 4 " 2m m _ e r(m-3/2 + e) 



= 2^" ^ T(m) y C(2m-3 + 2£ ), (A.3) 

rf-?§- = (^)-^-^- ^ m - 3 / 2 + e W - 5 + 26) (A 4) 

f p( p 2)m " l 47r 2 T2 J 22^-5^ r(m) 7 + 2e )- ^-5) 
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In appendix C, we give an example of how to calculate the above sum-integrals. More 
specifically we need 
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(A.10) 
(A.ll) 



We also need some fermionic one-loop sum-integrals. They can be obtained from the 
bosonic ones by scaling arguments |)2j]. The relations are 
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The specific sum-integrals needed are 
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(A.12) 
(A.13) 
(A. 14) 
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The two-loop integrals we need are 
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(A.21) 
(A.22) 



The last two-loop sum-integral can be obtained from the second by a change of variables. 
Moreover, in appendix C, we demonstrate how to calculate them. 



The simplest three-loop sum diagrams are the bosonic, fermionic and mixed basketballs. 
They read 



1 



pqk P 2 Q 2 K 2 (P + K + Q) 2 



1 /T 2 N 2 

(4tt) 2 V 12 



6 , a 182 



-12 



CM) , , C'(-i) 



{pqk} P 2 Q 2 K 2 (P + K + Q) 2 



C(-3) 
tt) 2 V12^ 



C(-i) 



+ 0(e), 



(A.23) 



3 , a 173 
2^ +9ln 4^+120 



t 



(Air) 

-^ln2-3^ + 12^ 
5 3 C("3) +12 C(-1) 



pq{k } P 2 Q 2 K 2 (P + K + Qy 



-(-)'' 
tt) 2 V12^ 



3 

4~e 



»h.Ji- 

2 4irT 



+ 0(e),(A.24) 
179 



(47T) 

+ ^ln2 + ^^l-6^ 
+ 10 ln2+ 2 C (-3) 6 C(-1) 



40 

+ 0(e). (A.25) 



In appendix C, we outline the calculation of the fermionic basketball. Finally, there are 
some more difficult three-loop sum-integrals: 



t 



(QK) 



1 /T 2 x2 



{p}qk P 2 Q 2 K 2 (P + K) 2 (P + Qf 



(QK) 2 

p{qk } P 4 Q 2 K 2 (P + K) 2 (P + Qf 



(Q - Kf 
pqk P*Q 2 K 2 (P + K) 2 (P + Q) 2 



3 9, u 9 361 
h -ln-^— + -7 E H 

8e 4 4ttT 4 [£ 160 



(47T) 

_^ln2 + ^-^ 
10 2C(-3) 2 C(-l) 

'T 2 x2 



+ 0(eXA-26) 



Itt) 2 V 12^ 



5 5, a 1 23 
h - In -!— + -7 B H 

24e 4 4vrT 4 ' 24 



8 In 2 1C/( ~ 3) I ?C ' ( ~ i: 
5 6 C(-3) + 6 C(-l) 



11 



+ 66 In 



+ C(e),(A.27) 
73 



3(4tt) 2 V12 

_ 10 ^ + 6 4-^ 



4^T + Y + 12 ^ 



C(-3) C(-i) 



+ 0(e) 



(A.28) 
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Appendix B 

Integrals in the Effective Theory 



In this appendix we define the integrals we need in the calculations in the effective three 
dimensional Euclidean field theory. We employ dimensional regularization in 3 — 2e di- 
mensions to regularize infrared and ultraviolet divergences. In analogy with Appendix A, 
we define 

(B.l) 



_ /e>Vv f d3 ~ 2e P r , 



, P - 4vr ) J (2?r) 

Again /x coincides with the renormalization scale in the modified minimal subtraction 
renormalization scheme. 



In the effective theory we need the following one-loop integrals 

fi 8 



ln(p + m 



P p 2 + m 2 
1 

p (p 2 + m 2 ) 2 



nr 

6ir 
m 

47T 

1 



57rm 



"l + (21n-^ + ^) e + 0(e 

v 2m 3 ' 

l + (21n-^ + 2)e + 0(e 2 
v 2m ' 

l + (21n^)e + 0(e 5 



(B.2) 
(B.3) 
(B.4) 



All integrals are straightforward to evaluate in dimensional regularization. Details may 
be found in Ref. [p|. 



The specific two-loop integrals needed are 

1 



P q (p 2 + mf) (q 2 + ml) [(p + q) 2 + m|] 



L 



pq (p 2 + m 2 ) 2 (q 2 + m 2 )(p — q) 2 



(4vr) 2 
In 



1 1 



mi + m 2 + m 3 



(47r) 2 m 2 



- A + 0(e) 



0(e, 



(B.5) 
(B.6) 
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Ipq (p 2 + M 2 )(g 2 + M 2 )[(p + q- 



-4 In 



2(M + m) 



k) 2 

M 



in 



rl 

- + 6 

e 



+ 4 — In , r 

m M + m 



+ 0(e) 



(B.7) 



The first two of these integrals can be found in Refs. [36,73]. The integral in Eq. ( |B.7| ) has 
previously been calculated by Braaten and Nieto for m = M in Ref. [58|. In appendix C, 
we calculate it for the more general case m ^ M. 



Appendix C 

Some Sample Calculations 



This appendix is devoted to the explicit calculation of some sum-integrals in order to 
illustrate the methods invented by Arnold and Zhai We will give some examples of 



how one computes one, two and three-loop diagrams and choose to present the calcula- 
tions of sum-integrals that have not explicitly worked out in the literature. 

Let us start with the following one-loop sum-integral: 

Po /e 7 Vv ^v- f d3 ~ 2t P (2t™T) 4 



p (_p2)m V 4tt / (2vr) 3 - 2e [p 2 + (27mT) 2 ]™ 

/e 7B /x 2 y 4 r d 3 ~ 2t p 1 ^ 7 



IE ,,2 r r 73-2e i 

(27r) 3 " 2e [p 2 + (2ttT) 2 

In the second line we have changed variables, and in the last line we have employed the 
definitions of the T and (^-functions, as well as performing a standard one-integral using 



dimensional regularization [100]. 



The next sum-integral we consider comes from the fermionic setting sun diagram: 

(C.2) 



i _ j n,(p) 



f{PQ} P 2 Q 2 (P + Q) 2 T P P 2 ' 
Here, we have defined the fermionic self-energy 

n ' (f H«wpW- (C3) 

There are different ways to compute this sum-integral. One is the contour method. One 
rewrites the sum over p and q as contour integrals. One then finds terms independent, 
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linear and quadratic in the distribution functions, which separately can be treated using 
dimensional regularization [f43| . We shall use the method of Arnold and Zhai, who have 



developed a new machinery to computing difficult multi-loop sum-integrals [42 



The idea is as follows. The calculation of an n-loop sum-integral obviously requires 
the evaluation of n sums and n three-dimensional integrations, if we work in momentum 
space. In coordinate space it only requires one four dimensional integration. At first sight, 
it therefore seems that it would be simpler to compute the sum-integrals in coordinate 
space. This would certainly be the case if the expressions were finite, so that we could 
evaluate the expressions directly in four dimensions. However, the sum-integrals are UV- 
divergent, and we must subtract off these divergences that arise at T = 0. This is most 
easily carried out in momentum space. The remainder may then be evaluated in four 
dimensions, and this is then done using the Fourier transform of the momentum space 
propagator: 



-\Po\r 



4nr 



(C.4) 



First, we separate the fermionic self-energy into a T = piece and a finite temperature 
term by writing 

U f (P) = uf\P) + uf\P). (C.5) 
The temperature independent part is, using standard results from dimensional regular- 



ization 100 



n (o) 



By using Eq. 



2c 



d d q 1 _ 1 4vr/i 2 

(2n) d Q 2 (P + Q) 2 ~ (47r) 2 ^~p2" 

one finds 



- lE + 2 + 0(e) 



(C.6) 



nf (P) 

p2 



1 T 2 /l 



(4tt) 2 12 V e 



(- + 4 In 



AixT 



6 + 4 



C(-i) 



(C.7) 



We now need TIj (P), and in order to obtain it we shall compute 11/ (P) and subtract off 
its T = limit. The fermionic self-energy is 



U f (P) = TY,Jd 3 rA(q ,r)A(q +p ,ry p - r 

d 3 r 



T 
(4^f 



The sum over q Q is given by 



y) e '" ' ( 

{go} 



ko\r „-\po+qo\r 



£ +' 



-\lo\r p -\po+qo\r Jp-r 



(C.8) 



[cschr + |po|K 



-\po\r 



(C.9) 



3.0 Some Sample Calculations 



121 



Here, f = 2irrT and po = po/2irrT. This formula can obtained by splitting the sum into 
three parts depending on the sign of go and (go + Po), and then use known results for 
geometric series. Substituting Eq. ( |C.9|) into Eq. ( |C.8|) and letting T — > 0, we obtain 

r (o; 

7 



nf\P). Subtracting this from Tlf(P), we find 



U T f (P) = — ^ / ^e ip ' r [cschf - l/f]e-l p °l r + 0(e). (CIO) 
3 (An) 2 J r 2 

Although IIj (P) is finite for P — ► oo, nJ(P)/P 2 is logarithmically divergent, because 
the former goes like 1/P 2 in this limit: 

This behaviour can be inferred by the contour trick: One rewrites the sum over Matsubara 
frequencies as a contour integral and study its high P limit [55]. Moreover, the high 
momentum behaviour is also given by the small r behaviour of the integrand in Eq. ( |U.10|) , 
and using the series expansion of cschf one finds 

Ilf \P) - --_ / —fe**e-Mr. (C.12) 

We can now write the finite temperature part of the fermionic setting-sum diagram as 

T 2 r 1 r r/ 3 r r 1 

(£54 ^ / ^ [CSChf - 1/f + r -/6]e'-eH«l' + 2^ {q) — r (C.13) 

Consider the first term above, which we denote by J. Integration over p simply gives the 
propagator in coordinate space, while the summation over p yields a factor cothf, since 

]T = e- |9o|r e- |po+(?o|r = [cothf + |po|]e- |po|r . (C.14) 

Furthermore, integration over the sphere gives the usual 4n, and we find 

T 2 r f dr i 

/ — [cschr — 1/f + r/6] cothf. (C.15) 
J r i 



(Any 

The integral that has been obtained are convergent. One can then calculate it numerically. 
However, it is divergent term by term, but Arnold and Zhai have developed a clever way 
to compute them analytically, using methods similar to dimensional regularization [^ . 
Below we shall discuss the derivation of them. 



The integral can be expressed in terms of T-functions and (^-functions: 



T 2 



(4vr)< 



(2 - T 



2- z - x )Y(z)i(z) + -2- z+1 T{z + 2)C(z + 2) 
o 



-1. (C.16) 
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We proceed by expanding the T and ^-functions around the pole z = — 1. This produces: 

(C.17) 



1 T 2 r C(-V 

2 lE -4 - 21n2 - 2- s v 



(4tt) 2 12 



C(-i) 



Finally, we have the second term in Eq. ( |C13| ), which is easily evaluated using appendix 
A: 



Tpw.\ PK 



1 T2 ^ 41n 44 + 2 ^ + 2 - 21n2 + 2 fFiy 



(C.18) 



fp {Q} P*Q 2 (47r )2 12 
Adding up the different pieces, we conclude that the fermionic setting sun graph vanishes: 



1 



f{p Q} PH) Q) 1 



0. 



(C.19) 



Let us now move on and discuss the evaluation of the mixed boson-fermion basketball 
diagram. We shall be slightly more sketchy this time. The sum-integral reads 



j = i 



i 



jpq{k} P 2 Q 2 K 2 (P + Q + K) 2 
By changing variables this may be rewritten as 

f p U f (P)U b (P). 

Here, IT^P) is the bosonic self-energy 

1 



(C.20) 



(C.21) 



r Q Q 2 {P + Q) 2 ' 

In complete analogy with the fermionic case, we write 

n 6 (P) = ni 0) (P) + nf ) (P). 

Note that T\f\P) = uf ] (P). The mixed basketball then reads 

j = ^nf (P)ni 0) (P) + ^nf ) (P)ni 0) (P) + ^nf (P)nf } (P) 
+^ p nf ) (P)nf ) (P). 

By using Eq (|A.3|) once more we derive the result 



(C.22) 



(C.23) 



(C.24) 



2_ , 12 24 12 C'(-3) 
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The next two terms are calculated using similar methods |42 |. In order not to get 



overwhelmed by calculational details, we simply state the result [(42 

-T 2 \2 



(0), 



-(-Y 



(4tt) 2 V12 



1 3 , a 
1 In- A 
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20e 10 4?rT 120 



n^(P)ni T) (P) 



_37 ln2 _ 37^31 + A Cp) 

io ioc(-3) C(-i) 
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(C.26) 
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4 24, a 103 
1 In H 

5e 5 4vrT 15 



+ 



+ 4- 



5 C(-3) C(-l) 
The bosonic self-energy may be written as 

T 



MP) 



A(qo,r)A(q + Po,ryr- r . 



(C.27) 



(C.28) 



As in the previous calculation, we find (P) by subtracting off its T = part 

(4^2 / ^"[oathf - iMe-'^. (C.29) 



n 



(T) 



The high momentum behaviour of the bosonic self-energy is obtained in analogy with the 
Eqs. ( |C. 11| ) and ( |C.12| ) in the fermionic case: 



nf } (P) 
nf } (P) 



1 T 



d 3 7 



3 (4tt) 2 

-H - 

P 2 T Q Q 2 



-fe ip ' r e~' Po ' r . 



(C.30) 
(C.31) 



The last term in Eq. (|C.24|) can then be rewritten as 

T 4 f df r _ _ _ _ f 2 

^ - J — (cschf — l/f)(cothf — l/f)(cothf — 1) + — (cothf — 1 



+ 4 



PQ{i<] P 4 Q 2 K 2 ' 
(C.32) 

Here we have performed the sum over Matsubara frequencies using Eq. ( U.14Q . This 
first integral above is again finite, but infinite term by term. Using the regularization 
techniques of Arnold and Zhai the first term above yields 



(4tt) 2 V12 
while the second reads 



1 /T 2 x2r 169 24, C'(-3) C'(-l) 

+ — In 2 + 2-4 — r-64 — r 

1 30 5 C(-3) C(-l) 



(47T) 



( — ) 

7r) 2 Vl2^ 



T 2 x2 r 2 A C'(-l) 

_ + 121n _ + 47E + 8-41n2 + 8^- I i 



(C.33) 



(C.34) 
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Adding all the terms, we obtain Eq. ( [A.25Q . 



We would also like to mention that it has been noted that one can obtain this diagram 



from the bosonic and fermionic basketball diagrams by scaling arguments ||43|| . One finds 

fp Q {K} P 2 Q 2 K 2 (P + K + Q) 2 = ~6 [1 ~ 2 3 ^kq P 2 Q 2 K 2 (P + K + Q) 2 

~6T{pqk} P 2 Q 2 K 2 (P + K + Qf ( °' 35) 

We have checked that our results agree. 

Finally, we must consider the calculation of the divergent integrals that appeared above. 
The first needed is 

con 

drr z . (C.36) 



o 



Depending on the value of z, the contribution to the integral from one of the limits van- 
ishes (and the other blows up). Thus, if one analytically continue z independently to 
regulate the behaviour at r = and r = oo, the integral in Eq. ( |C.36| ) vanishes. 



We also have the result 

drr z e- ar = a- 1 - z T(l + z), (C.37) 



oo 

ar „ — 1— 2i 







for the values of z for which the integral is well defined. The integral is then defined 
for all values of z by analytic continuation. This makes it possible to attack integrals of 
hyperbolic functions times powers of r: 



poo poo r ^ 

/ drr z cschr = r z 2e~ r 
Jo Jo - „_ n 



oo 

e ~2nr 



n=0 

oo 



2-r(i+*)i; 1 



= (2-2~ z )V(l + z)C(l + z). (C.38) 
In the last line, we have used the functional relation between the Riemann Zeta-function 



and the Hurwitz Zeta function ITS 



In a similar fashion one can obtain 



POO 

/ drr z cothr = 2~ z T{z + 1)C0+ !)• (C.40) 
Jo 
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More complicated integrals, such as products between cschf and coth f can be computed 
from the above formulas using integration by parts. We then have all regulated integrals 
needed to evaluate the expressions above. 



Let us now turn to the three-dimensional integrals. All, except for one integral in the 
effective have been computed by Braaten and Nieto in Ref. In order to illustrate 

their methods, we shall calculate this integral, which is 



pq (p 2 + M 2 )(q 2 + M 2 )[(p + q + k) 2 + m 2 



k=im 



(C.41) 
It can best be 



The integral has been computed in the less general case m = M 
computed by going to coordinate space. The Fourier transform of the propagator is 

1 



V m (R) 



p 2 + m 2 

It can be expressed in terms of a modified Bessel function 



V m (R) = ( 



Att > (2n)V 2 -^R 



I /m\ i/2-« 

— (17) K 1/2 _ e {mR). 



In three dimensions (e = 0) this is the Yukawa potential: 

p —mR 

Vm(R) 



AtcR 



(C.42) 



(C.43) 



(C.44) 



For small R it can be written as a sum of two Laurent series in R 2 . One of these is singular 
beginning with an R~ 1+2e term and the other is regular which begins with an R° term: 



V m {R) 



^/A.ra-e) 1 



rfi) 4vr 



R 



-l+2e 



m 2 R 2 



An 



1 + 



2(1 + 2e) 
m 2 R 2 



2(3 - 2e) 



+ 0(m 4 i? 4 ) 
+ 0{m i R A 



(C.45) 
(C.46) 



The integral can be written 



1 



/ ^V 2 M {R)V m {R). (C.47) 

JR 



l P q (p 2 + M 2 )(q 2 + M 2 ) [(p + q + k) 2 + m 2 ] 

The radial integration is now split into two regions, < R < r and r < R < oo. The 
ultraviolet divergences arise from the region R — > 0. This implies that we can set e = 
in the region where r < R < oo. Hence, one can write the integral as 

-e^Vx-e^Tr) 3 / 2 rr 



e^V^(R)V m (R) 



+■ 



2k 

An r°° 
k 



dRR^J l/2 _ t {kR) Vl (R) V m (R) 



\fk Jo 
dRR sm(kR)V^(R)V m (R) 



(C.48) 
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Here, J v (x) is an ordinary Bessel function. The Bessel function has the following expansion 
for small R: 

'hR) 1/2 ~\l + 0{k 2 R 2 )]. 



Jl/2-e{kR) 



r(f-e)V2 



(C.49) 



Using this expansion and the small R expansion of the potential, the first integral is, after 
dropping terms that vanish in the limit r — > 



V 2k > 



(2tt 



,3/2 



Vk 



dRR 3 / 2 -tJ 1/2 ^(kR)V^(R)V m (R) 
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•57T 



— h 4 In /ir 



+2 + 4 7i; +O(eXC50) 



The second integral can be found in e.g Ref. [15| and equals 



(2M+m+ik)T[-l,(2M+m+ik)r]-(2M+m-ik)T[(2M+rn-ik)r] . (C.51) 



2fc(47r) 2 

Evaluating this at k = im yields 

— / dIffiam(kR)Vu(R)V m (R) 
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m M + m 



- 7s + 1 



-ln[2(M + m)r] + 0(e) 
where we have used the series expansion of the incomplete gamma function 



(C.52) 



r[-l,x] = - + 7B- 1 + lna; + 0(x 2 ), 



(C.53) 



and dropped terms that vanish as r — > 0. Collecting our results we obtain Eq. ( |B.7| ). The 
logarithms of r cancel and our result reduces to the one found in Ref. in the case 
m = M, as it should. 
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